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ABSTRACT 

High-redshift submillimetre galaxies (SMGs) are some of the most rapidly star-forming galax- 
ies in the Universe. Historically, galaxy formation models have had difficulty explaining the 
observed number counts of SMGs. We combine a semi-empirical model with 3-D hydrody- 
Q^' namical simulations and 3-D dust radiative transfer to predict the number counts of unlensed 

SMGs. Because the stellar mass functions, gas and dust masses, and sizes of our galaxies 
£^ ■ are constrained to match observations, we can isolate uncertainties related to the dynamical 

^ | evolution of galaxy mergers and the dust radiative transfer. The number counts and redshift 

Ci ■ distributions predicted by our model agree well with observations. Isolated disc galaxies dom- 

inate the faint (5i.i < 1 mJy, or Ssso < 2 mJy) population. The brighter sources are a mix 
of merger-induced starbursts and galaxy-pair SMGs; the latter subpopulation accounts for 
~ 30 — 50 per cent of all SMGs at all Si.i > 0.5 mJy (S'sso 1 mJy). The mean redshifts are 
£^ " ~ 3.0 — 3.5, depending on the flux cut, and the brightest sources tend to be at higher redshifts. 

Because the galaxy-pair SMGs will be resolved into multiple fainter sources by ALMA, the 
bright ALMA counts should be as much as 2 times less than those observed using single-dish 
fS| ' telescopes. The agreement between our model, which uses a Kroupa IMF, and observations 

suggests that the IMF in high-redshifts starbursts need not be top-heavy; if the IMF were 
■ top-heavy, our model would over-predict the number counts. We conclude that the difficulty 

| some models have reproducing the observed SMG counts is likely indicative of more general 

problems - such as an under-prediction of the abundance of massive galaxies or a star forma- 
tion rate-stellar mass relation normalisation lower than that observed - rather than a problem 
specific to the SMG population. 
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1 INTRODUCTION 



Submillimetre galaxies (SMGs; Smail et ail 1 19971: Barger et al 
1 19981 : iHughes et alJl 19981 : lEales et al] 1 19991 : see iBlain et al J 12002 
for a review) are amongst the most luminous, rapidly star-forming 
galaxies known, with luminosities in excess of 10 12 L and 
star f ormation rates (SFR) of order ~ 10 2 - 10 3 M yr" 1 
(e.g., iKovacs et al] 120061 : ICoppin et al] |2008t IMichalowski et al] 
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|201CL 120121 : iMagnelli et all |201CL 1201 2l : IChapman et ail l2010h . 
They h ave stellar masses o f ~ 10 11 Mr, although recent esti- 
mates i lHainline et al] I20T1I : IMichalowski etal J |20 id 120121) dif- 
fer by a factor of ~ 6, and typical gas fract i ons o f ~ 40 



per cent llGreve et al] 120051: 



iNaravanan. Bothwell. & Pavel 



Tacconi et al] 120061, 120081 : but cf. 
2012al) . 



The most luminous local galaxies, ultra-luminous infrared 
galaxies (ULIRGs, defined by L m > 10 12 L B ), are almost ex- 
clusively late-stage major mergers (e.g., lLonsdale et al]|2006t) be- 
cause the strong tidal torques exerted by the galaxies upon one 
another when they are near coalescence cause signific ant gas in- 
flows and, consequently, bursts of star formation (e.g.. iHernquistl 
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1 19891 : lBarnes&Hemguisllll99lL 1 19961 : iMihos & Hemguistlll996t) . 
Thus, it is natural to suppose that SMGs, which are the most 
luminous, highly star-forming galaxies at high redshift, are also 
late-stage major mergers undergoing starburst s. There is signifi- 
cant observa t ional support for this picture (e.g. Jlvisqn et alj|2002t 
200 1 1201(1 IChapman et al l |2003t iNeri et alj [20031: Smail et alj 



2004l;ISwinbank et alfcOofereve et 



: al J2005I: 



Tacconi et al. 2006, 



Younger et al. 


2008 


, 201oT Iono et al. 20091: lEngel et al 201C 


Bothwell et alj 


201C 


,l2012l;lRiechers et alj|20 11 alibi Magnelli et al 



ies of t he required masse s to account for the observed SMG abun- 
dances dDave et al.l2010l) . Consequently, explaining the abundance 
of SMGs has proven to be a challenge for galaxy formation models. 

Much observational effort has been invested to determine 
the number counts and redshift distribution of SMGs (e.g., 



Chapman et alj |2005|; ICocoin et alj 


20061: Knudsen et alj 


2008: 


ChaoinetalJ 20091: IweiBetalJ 20C 


9l: lAustermann et alj 


2005 


20ld; Scott et al.Eoiof Zemcov et al 


1 20 id: Aretxaga et al. 


2011; 


Banerii et alj 201 1 ; lHatsukade et alj 


EOllt IWardlowetal 


2011; 


Roseboom et alj|20i2t lYun et all 1201 


2|) because this information 



is key to relate the SMG population to their descendants and 
to understand SMGs in the context of hierarchical galaxy for- 
mation models. Various authors have attempted to explain the 
observed abundance of SMGs using phenomenological mod- 
els (e.g.. |Pearson&Rowan-RobinsorJll996t iBlain et alj Il999bl: 



Devriendt & Guiderdonill200d:lLagache et al.ll2003l ; lNegrello et alj 



20071 : iBefhermin et all |2012|), semi-ana l ytic m o dels (SAMs; 
e.g., iGuiderdoni et alj Il998t IBlain et al] 1999tJ ; iGranato et aH 



2000; Kaviani et al. 120031; IG ranato e t alj|2004l; iBaugh et alj 200f 
Fontanot et alj|20 07; Fontanot & Mo nacdl20ld:lLacev et alj|200fc 



2011 



2001 



20ld: ISwinbank et alj|2008t ILq Faro et alj 120091: iGonzalez et al 



) , and cosmological hy drodynamical simulations ( Fardal et al 
Dekel et all2009tlDave et alj|20ld : lshimizu et alj|2012h . 



Granato et alj ( 12000!) presented one of the first SAMs to self- 
consistently calcul ate dust absorptio n and emission by coupling the 
GALFORM SA M dCole et alj|200d) with the GRASIL spectropho- 
tometric code dSilva et al.lll998h . This was a significant advance 
over previous work, which effectively treated the dust temperature 
as a free parameter. Self-consistently computing dust temperatures 
made matching the submm coun ts significantly more difficult: the 
submm counts predicted by the IGranato et alj model were a fac- 
tor of ~ 20 — 30 less than those observed jBaughetal.ll2005l ; 
ISwinbank et"ai]|2008l). 

The work of B augh et a l] J2005I. hereafter B05) has attracted 
significant attention to the field because of its claim that a flat IMF 
is necessary to reproduce the properties of the SMG population, 
which we wil l discu ss in detail here. B05 set out to modify the 
IGranato et al. I d2000h model so that it would reproduce the prop- 
erties of both 2 ~ 2 SMGs and Lyman-break galaxies (LBGs) 
while also matching the observed 2 = optical and IR luminos- 
ity functions. Ad opting a flat IMlf^] in starbursts ra ther than the 
iKennicuttl l ll983h IMF used in IGranato etai] l l200d) was the key 
change that enabled B05 to match the observed SMG counts and 
redshift distribution while still reproducing the local A'-band lu- 
minosity function. A more top-heavy IMF results in both more 



1 Specifically, the IMF they use i s dn / dlog M = constant for the mass 
range 0.15 < M < 125 Mq, The lKroupall200ll) IMF has dn/dlos M oc 
M -1 - 3 for M > 1 Mq, so the difference between the B05 IMF and that 
observed locally is considerable. 



luminosity emitted and more dust produced per unit SFR; conse- 
quently, the s ubmm flux per unit S FR is increased significantly 
(see B05 and lHavward et all 1201 lal . hereafter Hll, for details). 
The B05 modifications increased the Ssso per unit SFR for star- 
bursts by a factor of ~ 5 (G.-L. Granato, private communication), 
which caused starbursts to account for a factor of ~ 10 3 t imes 
more sources at Ssso = 3 mly than in IGranato et al] J2000h . As 
a result, in the B05 model, ongoing starbursts dominate the counts 
for 0.1 < Ssso i$ 30 mjy. Interestingly, thes e starbursts are 
trigge r ed predominantl y by m inor mergers (B05; [Gonzalez et alj 
l201ll) . ISwinbank etai] d2008h present a detailed comparison of 
the properties of SMGs in the B05 model with those of observed 
SMGs. The far-IR SEDs, ve locity dispersions, and halo masses 
(see also lAlmeida et al]|201 ll) are in good agreement; however, re- 
cent observations suggest that the typical redshift of SMGs may be 
higher than predicted by the B05 model and, contrary t o the B05 
predi c tion, brighter SMGs tend to be at higher redshifts dYun et al] 
120121 ; ISmolcic et alj 12012]) . Furthermore, the rest-frame A'-band 
fluxes of the B05 SMGs are a factor of ~ 10 lower than observed; 
the most plausible explan ation is that the mass es of the SMGs in 
the B05 SAM are too low dSwinbank et al.l2008h . but the top-heavy 
IMF in starbursts used by B05 makes a direct comparison of masses 
difficult. These disagreements are reasons it is worthwhile to ex- 
plor e alternative SMG mo dels. 

IGranato et alj d2004l) presented an alternate model, based on 
spheroid formation via monolithic collapse, that predicts submm 
counts in good agreement with those observed and reproduces the 
evolution of the A-band luminosity function. However, the typical 
redshift the y predict for SMGs is lower than rec ent observational 
constraints dYun et alj [20121 : ISmolcic et alj|2012l) , and this model 
does not include halo or g alaxy m ergers. 



The iFontanot etalT l l2007l) model predicts SMG number 
counts in reasonable agreement with those observed using a stan- 
dard IMF; they argue that the crucial difference between their 
model and that of B05 is the cooling model used (see also 

1 Viola etalj|200i and |Pe Lucia et alj|20ich . However, their SMG 
redshift distribution peaks at a lower redshift than the redshift 
distribution derived from recent ob servations ( Yun et al] [20121 : 
ISmolcic et al.l2012l) . Furthermore, the lFontanot et alj J2007) model 
produces an overabundance of bright galaxies at z < 1. How- 
ever, this problem has been significantl y reduced in the latest ver- 
sion of the model dLo Faro et al.ll2009T) . which provides a signifi- 
cantly better fit to th e galaxy stellar mass function at low redshift 
dFontanot et al.ll2009T) . In the revised model, the submm counts are 
reduced by ~ 0.5 dex, primarily because of the change in the IMF 
from Salpeter to Chabrier, but the redshift distribution is unaffected. 
Thus, the submm counts for the new model are consistent with the 
data for Ssso S; 3 mjy, but they are slightly less th an the observed 
counts at higher fluxes dFontanot & Monacal 20 id) . No fine-tuning 
of the dust parameters has been performed for the new model. 

A compelling reason to model the SMG population in an alter- 
native manner is to test whether a top-heavy IMF is required to ex- 
plain the observed SMG counts. Matching the submm counts is the 
primary reason B05 needed to adop t a flat IMF in starbursts^ Using 
the same model. lLacev et all d2008l) show that the flat IMF is neces- 
sary to repr oduce the evolution of the mid-IR luminosity function. 
Others (e.g jGuiderdoni et alJll998tlBlain et aHll999d ; |Pave et alj 

2 Recent observations suggest the number counts are a s much as a factor 
of 2 lower than those used by B05. Thus, if lBaugh et alj were to attempt to 
match the revised counts, the required IMF variation would be more modest. 
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have also suggested that the IMF may be top-heavy in SMGs, 
but they do not necessarily require variation as extreme as that as- 
sumed in B05. However, the use of a flat IMF in starbursts remains 
controversial: though there are some theoretic al reasons to believe 
the IMF is more top-heavy in starbursts (e.g., Larsonl 199S 
Elmegreen & Shadmehril 120031 ; lElmegreenl 120041 : iHopkim 



2005 



2012 



Narayanan & DaveT 20121) . there is to date no clear evidence for 



strong , systematic IMF variation in any environment l lBastian et alj 
|2010| and references therein). Furthermore, in local massive el- 
lipticals, the probable desce ndants of SMGs, the IMF may ac- 
tually be bottom-heavy (e.g., van Dokkum & Conroy||2010l . l201ll ; 



IConrov & van Dokkuml2012| - |Hopkinsl20'l2Tr Finally, the large pa- 
rameter space of SAMs can yield multiple quali tatively distinct so- 
lutions that satisfy al l observational constraints dBower et alj20ld : 
iLu et al ]|201ll . l2012l) . so it is possible that a top-heavy IMF in star- 
bursts is not required to match the observed submm counts even 
though it enables B05 to match the submm counts0Thus, it is use- 
ful to explore other methods to predict the submm counts and to 
determine whether a match can be achieved without using a top- 
heavy IMF. 

Another reason to model the SMG population is to in- 
vestigate whether, like local ULIRGs, they are predominantly 
merger-induced starbursts. Some observational evidence suggests 
that some SMGs may be early-stage mergers in which the discs 
have not yet coalesced and are likely not undergoing starbursts 
(e.g. iTacconi et alj|2006l. 120081 ; lEngel et alj[201ol : iBothwell et alj 
l2O10t iRiechers et alj 1201 1 alibi) , and massive isola ted disc galax- 
ies may also contribute t o the population ( e.g. , iBothwell et al 



2010 j: (Carilli et alj 120101: iRicciardelli et alj 120101 : iTargett et al . 
201 ll l2012h In Hll and HI 2, we suggested that the inefficient 



scaling of (sub)mm flux with SFR in starbursts results in an SMG 
population that is heterogeneous: major mergers contribute both as 
coalescence-induced starbursts and during the pre-coalescence in- 
fall stage, when the merging discs are blended into one (sub)mm 
source because of the large (~ 15", or ~ 130 kpc at z ~ 2 — 3) 
beams of the single-dish (sub)mm telescopes used to perform large 
SMG surveys. We refer to the latter subpopulation as 'galaxy-pair 
SMGs'. Similarly, compact groups may be blended into one source 
and can thus also contribute to the population. The most mas- 
sive, highly star-forming isolated discs may also contribute (Hll). 
Finally, it has been observationally demonstrated that there is a 
contrib ution from physic ally unrelated galaxies blended into one 
source dWang et alJioTTh . It is becoming increasingly clear that the 
SMG population is a mix of various classes of sources; if one sub- 
population does not dominate the population, physically interpret- 
ing observations of SMGs is significantly more complicated than 
previously assumed. 

In previous work, we demonstrated that major merg- 
ers c an reproduce the obse rved 850-|i.m fluxes and typical 
SED dNaravanan etall l2010bl): CO spatial exte nts, line-widths, 
and excitation ladders 1 Narayanan et alj [20091); s tellar masses 
dNaravanan et alj l2010bl : Hll; iMichaiowski et alj |2012|) : L m - 
effective dust temperature relation, IR excess, and star forma- 
tion efficiency dHavward et alj|2012l . hereafter HI 2) observed for 
SMGs. In this work, we present a novel method to predict the 
(sub)mm counts from mergers and quiescently star-forming disc 
galaxies. We utilise a combination of 3-D hydrodynamical simula- 
tions, on which we perform radiative transfer in post-processing to 

3 However, whether a unique solution, if any solution at all, can be found 
when all possible observational constraints are included is an open question. 



calculate the UV-to-mm SEDs, and a semi-empirical model (SEM) 
of galaxy formation - both of which have been extensively vali- 
dated in previous work - to predict the number counts and red- 
shift distribution of SMGs in our model. We address four primary 
questions: 1. Can our model reproduce the observed SMG number 
counts and redshift distribution? 2. What are the relative contribu- 
tions of merger-induced starbursts, galaxy pairs, and isolated discs 
to the SMG population? 3. How will the number counts and red- 
shift distribution of ALMA-detected SMGs differ from those de- 
termined using single-dish surveys? 4. Does the SMG population 
provide evidence for a top-heavy IMF in high-redshift starbursts? 

The remainder of this paper is organised as follows: In Section 
[3] we present the details of the simulations we use to determine the 
time evolution of galaxy mergers and to translate physical proper- 
ties of model galaxies into observed-frame (sub)mm flux densities. 
In Section [4] we discuss how we combine the simulations with a 
SEM to predict the (sub)mm counts for merger-induced starburst 
SMGs (Section [4. It and isolated disc and galaxy-pair SMGs (Sec- 
tion !4.2b . In Section[5] we present the predicted counts and redshift 
distribution of our model SMGs and the relative contribution of 
each subpopulation. We discuss implications for the IMF, compare 
to previous work, and highlight some uncertainties in and limita- 
tions of our model in Section[6] and we conclude in Section|7] 



2 SUMMARY OF THE MODEL 

Predicting SMG counts requires three main ingredients: 1. Because 
SFR and dust mass are the most important properties for predicting 
the (sub)mm flux of a galaxy (Hll), one must model the time evo- 
lution of those properties for individual discs and mergers. 2. The 
physical properties of the model galaxies must be used to deter- 
mine the observed-frame (sub)mm flux density of those galaxies. 
3. One must put the model galaxies in a cosmological context. Ide- 
ally, one could combine a cosmological hydrodynamical simulation 
with dust radiative transfer to self-consistently predict the (sub)mm 
counts. However, this is currently infeasible because the resolu- 
tion required for reliable radiative transfer calculations cannot be 
achieved for a cosmological simulatio n large enough to contain a 
significant number of SMGs (see, e.g jDave et alJ|201(il) FI 

Here, we develop a novel method to predict the number counts 
and redshift distribution of high-z SMGs while still resolving the 
dusty ISM on scales of ~ 200 pc. We predict (sub)mm c ounts using 
a combination of a simple SEM dHopkins et alj[2008al l3) and ide- 
alised high-resolution simulations of galaxy mergers. The method 
we use for each of the three model ingredients depends on the 
subpopulation being modelled. The physical properties of the iso- 
lated disc galaxies and early-stage mergers are determined using 
the SEM. For the late-stage mergers, hydrodynamical simulations 
are used because of the complexity of modelling a merger's evo- 
lution. Dust radiative transfer is performed on the hydrodynami- 
cal simulations to translate the physical properties into observed 
(sub)mm flux density. For the isolated discs and early-stage merg- 
ers, fitting functions derived from the simulations are used, whereas 



4 Recently, IShimizu et all fe012l) predicted SMG number counts using a 
cosmological simulation with a self-consistent model to calculate the far- 
IR emission. However, their model assumes a single dust temperature and 
neglects dust self-absorption, so the submm fluxes predicted by their model 
may be significantly (~ 0.3 — 0.5 dex) greater those calculated using full 
3-D dust radiative transfer (Hll). An investigation of the effects of this 
uncertainty on the predicted counts is underway. 
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Table 1. Summary of methods 



Ingredient Isolated discs Early-stage mergers Merger-induced starbursts 

Physical properties semi-empirical semi-empirical simulations 

(Sub)mm flux density HI 1 relations HI 1 relations simulations 

Cosmological context observed SMF merger rates from SEM merger rates from SEM 

+ duty cycle from sims + duty cycle from sims 



for the late-stage mergers, the (sub)mm light curves are taken di- 
rectly from the simulations. Finally, the isolated galaxies are put 
in a cosmological context using an observed stellar mass function 
(SMF). For the mergers, merger rates from the SEM and duty cy- 
cles from the simulations are used. The methods are summarised 
in Table [T] and each component of the model is discussed in detail 
below. 

We emphasize that we do not attempt to model the SMG pop- 
ulation in an ab initio manner. Instead, we construct our model so 
that the SMF, gas fractions, and metallicities are consistent with 
observations. This will enable us to test whether, given a demo- 
graphically accurate galaxy population, we are able to reproduce 
the SMG counts and redshift distribution. If we are not able to re- 
produce the counts and redshift distribution, then our simulations or 
radiative transfer calculations must be incomplete. If we can repro- 
duce the counts and redshift distribution, then it is possible that the 
failure of some SAMs and cosmological simulations to reproduce 
the SMG counts may be indicative of a more general problem with 
those models (e.g., a general under-prediction of the abundances of 
massive galaxies) rather than a problem specific to the SMG popu- 
lation. 

In the next two sections, we describe our model in detail. 
Readers whom are uninterested in the details of the methodology 
may wish to skip to Section[5] 



3 SIMULATION METHODOLOGY 
3.1 Hydrodynamical simulations 

We have performed a suite of sim ulations of isolated and merg- 
ing d isc galaxies wit h GADGET-2 dSpringel et alj|200lt ISpringell 
l2005h . a TreeSPH faernquist & Katzl ll989l) ~~code that com- 
putes gravitational i nteractions via a hierarchical tree method 
dB antes & Hull Il986l) and gas dy n amics via smoothed-particle 
hydrodynamic s (SPH; ILucvI 1 19771 : iGingold & Monaghanl 1 19771 ; 
ISpringell |2010|) FI I t explicitly conserves both e nergy and entropy 
when appropriate dSpringel & Hernguistll2002T) . Beyond the core 



5 Recently, some a u thors lAgeitz et al.1 l200l ISpringe! |2010| ; 
iBauer & Springe! 120121 ; ISiiacki et al.l |2012|) have noted several signif- 
icant flaws inherent in the traditional formulation of SPH, including the 
artificial suppression of fluid instabilities, artificial damping of turbulent 
eddies in the subsonic regime, and a lack of efficient gas stripping of 
infalling structures. Consequently, the results of cosmological simulations 
performed using GADGET-2 can differ significantly from those perfo rmed 
with the more accurate moving-mesh code AREPO jspringel 201 (j) even 
when the physics included in the codes i s identical fVoge lsberger et al. 
1201 it iKeres et aljlioil | Torrev et alj|201 ll) . Fortunately, a comparison of 
idealised merger simulations run with GADGET-2 and AREPO suggests 
that these issues do not significantly alter the global properties (e.g., star 
formation histories) of the mergers (Hayward et al., in preparation), so the 
results presented here should be robust to these numerical issues. 



gravitational and gas physics, the ver sion of GADGET - 2 we use in- 
cludes radiative heating and cooling dKatz et al]|l996t) . Star forma- 
tion is implemented using a volume-density-dependen t Kennicutt- 
Schmidt (KS) law l lSchmidtl 1 19591 ; iKennicuttl Il998l) , p S FR oc 
Pgas! with a low-density cutoff. We use N = 1.5, which re- 
produces the global KS law and is consistent with observa- 
tions of high-redshift disc ga laxies jKrumholz & Thompson |2007| ; 
iNaravanan et alJl2008Ll201 ll ; but see lNaravanan et alj|2012bf) . 

Furthermore, our simulations include a two-phase 
sub-resolution model for the interstellar medium (ISM; 
ISpringel & Hernquistl 120031) in which cold dense clouds are 
in pressure equilibrium with a diffuse hot medium. The division of 
mass, energy, and entropy between the two phases is affected by 
star formation, radiative heating and cooling, and supernova feed- 
back, which heats the diffuse phase and evaporates the cold clouds 
dCox et al.ll2006bl) . Metal enrichment is calculated assumi ng each 
particl e behaves as a closed box the yield appropriat e for a Kroupal 



2001 ) IMF. The simulations also include the ISpringel et alj 



2005) model for feedback from active galactic nuclei (AGN), 



in which black hole (BH) sink particles, initialised with mass 
lO 5 ^ -1 Mo, undergo Eddington-limited Bondi-Hoyle accretio n 



faovle & Lvttletonl 1 19391 ; iBondi & Hoylel 1 1944 iBondil Il952l) . 
They deposit 5 per cent of their luminosity (L = O.lrac 2 , where 
rh is the mass accretion rate and c is the speed of light) to the 
surrounding ISM. This choice is mad e so that the normalisa tion of 
the Mbh — a relation is recovered jPi Matteo et al.1120051) . Note 
that our results do not depend crucially on the implementation of 
BH accretion and feedback for two reasons: 1. the AGN typically 
do not dominate (but can still contribute significantly to) the 
luminosity of our model SMGs because the SEDs during the phase 
of strong A GN activity tend to b e hotter than during the starburst 
phase (e.g., lYounger et a"i]|2009l Snyder et al., in preparation), so 
the mergers are typically not SMGs during the AGN-dominated 
phase. 2. Even in the absence of AGN feedback, the SFR decreases 
sharply after the starburst simply because the majority of the cold 
gas is consumed in the starburst. 



m disc g a 

iHernquisl (l990) profile and an exponential gas and stellar disc 
in which gas initially accounts for 80 per cent of the total bary- 
onic mass. At merger coalescence, the baryonic gas fractions are 
typically 20-30 per cent , which is consistent with the estimates of 
INaravanan et alj ( l2012"3) . The mass of the baryonic component is 4 
per cent of the total. The galaxies are scaled to z = 3 following the 
method described in lRobertson et alj d2006t) . Dark matter particles 
have gravitational softening lengths of 200/i _1 pc, whereas gas and 
star particles have lOO/i -1 pc. We use 6 x 10 4 dark matter, 4 x 10 4 
stellar, 4 x 10 4 gas, and 1 BH particle per disc galaxy. The detailed 
properties of the progenitor galaxies are given in Table [2] Note 
that we have chosen galaxy masses such that most of the merg- 
ers, based upon our simulations, will contribute to the bright SMG 
population (i.e., at some time during the simulation they have ob- 
served 850-^m flux density Ssso > 3 mly). More massive galaxies 
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Table 2. Progenitor disc galaxy properties 





V 20 o fl 


A/200* 


c c 




-^gas,init 


R d f 


Name 


(kms^ 1 ) 


(lO 12 /^ 1 M ) 




(lO 10 ^- 1 Mq) 


(lO 10 ^- 1 M ) 


{hT 1 kpc) 


b6 


500 


6.4 


2.3 


5.3 


22 


3.3 


b5.5 


400 


3.3 


2.5 


2.7 


11 


2.6 


b5 


320 


1.7 


2.8 


1.4 


5.6 


2.0 


b4 


260 


0.60 


3.2 


0.49 


2.0 


1.7 



" Virial velocity. * Virial mass. c Halo concentration. d Initial stellar mass. e Initial gas mass. ' Disc 



scale length. 



Table 3. Merger parameters 



Name 


A*' 


R b 

-ii-pcri 

(h,- 1 kpc) 


■Rinit 

(h' 1 kpc) 


b6b6 


1 


6.7 


70 


b6b5.5 


0.52 


6.7 


70 


b6b5 


0.26 


6.7 


70 


b6b4 


0.09 


6.7 


70 


b5.5b5.5 


1 


5.3 


57 


b5b5 


1 


4.0 


44 



"Baryonic mass ratio Afb,secondary /-^primary- 
* Pericentric-passage distance. c Initial separation of 
the discs. 



will also contribute but are increasingly more rare, so our simula- 
tions should be representa tive of all but the brightest, rarest SMGs 
(Michatowski et al. IHoH. Note also that we have included some 
slightly lower mass mergers for completeness. 

We simulate each disc galaxy listed in Table [2] in isolation 
for l.bh^ 1 Gyr and use these isolated disc simulations as part of 
our simulation suite. Our suite also includes a number of simula- 
tions of major and minor galaxy mergers. For the merger simula- 
tions, two of the progenitor disc galaxies are placed on p arabolic 
orbits (which are motivated by c osmological simulations: iBensonl 
120051 : iKhochfar & Burked 120061) with initial separation i? init = 
5R200 /8 and pericentric-passa ge distance equal to twice the disc 
scale length, R pa ri = 2i?^ dRobertson et alj|2006t) . The evolution 
of the system is followed for 1.5/i _1 Gyr, which is sufficient time 
for the galaxies to coalescence and for significant star formation 
and AGN activity to cease. The details of the merger simulations 
are given in Table [3] For each combination of progenitor discs i n 
Table[3] we simulate a subset of the i-p orbits o flCoxetalJ ( f2006ah . 
Specifically, we use the i-p orbits for the major mergers (b6b6, 
b5.5b5.5, and b5b5) and the i and j orbits for the unequal-mass 
mergers (b6b5.5, b6b5, and b6b4) because the latter have shorter 
duty cycles and the variation in duty cycles amongst orbits is not a 
primary source of uncertainty. Consequently, we use a total of 34 
GADGET-2 simulations. 



3.2 Dust radiative transfer 

In post-processing, we use the 3-D Monte Carlo radiative trans- 
fer code SUNRISE to calculate the UV-to-mm SEDs of the 
simulated galaxies. We have previously si mulated galaxies with 
colours/SEDs co nsistent wit h local SINGS dKennicutt et all 20031 ; 
iDale et al.l |2007|) galaxies Jonsson. Groves. & Coxl l2010h ; lo- 
cal ULIRGs ( I Younger et all (20091) : massive, quiescent, compact 



SUNRISEis publicly available at jhttp :/ / code ■ google ■ com/p/ sunrise 



2 galaxies dWuvts et alj 120091 |2010|) : 24 nm-selected 
,alaxies 1 Narayanan et alj uOlOal) : K+A/post -starburst galaxie s 
Snyder et al.l 1201 lb : and extended UV discs dBush et alj boifjb . 
among other populations, so we are confident that SUNRISE can be 
used to model the high-z SMG population. As discussed above, pre- 
vious work has demonstrated that many properties of our s i mulated 
SMGs agree with observa tions dNaravanan et alj I2009L 1201 Obi : 
lHavward et alfeoi lall2012h , but we have yet to put our simulated 
SMGs in a cosmological context. We briefl y review the d e tails o f 
SUNRISE her e , but we re fer the reader to I Jonsson et alj d2006l) . 
Ijonsson et alj d201fjb . and I Jonsson & Primackl d201fjb for full de- 
tails of the Sunrise code. 

SUNRISE uses the output of the GADGET-2 simulations to 
specify the details of the radiative transfer problem to be solved, 
specifically the input radiation field and dust geometry. The star 
and BH particles from the GADGET-2 simulations are used as 
sources of emission. Star particles are assigned STARBURST99 
dLeitherer et~ai]| 1999b SEDs acc ording t o their ages and metallic- 
ities. We conservatively use the iKroupal d200lb IMF when calcu- 
lating the simple stellar population SED templates. Star particles 
present at the start of the GADGET-2 simulation are assigned ages 
by assuming that their stellar mass was formed at a constant rate 
equal to the star formation rate of the initial snapshot, but the re- 
sults are insensitive to this choice because we discard the early 
snapshots and the stars present at the start of the simulation ac- 
count for a small fraction of the luminosity at later ti mes. The initial 
gas a nd stellar metallicities are Z = 0.015 (~ Z® ; lAsplund et al] 
2009)- We have chosen this value so that the starbursts lie roughly 
on the observed mass-metallicity relation; however, the results are 
fairly robust to this choice because a factor of 2 change in dust 
mass changes the (sub)mm flux by only ~ 50 per cent because 
(sub)mm flux scales approximately as M^ 6 (Equations! 1 5 land[T6t . 
BH particles are assigned luminosity-depen dent templates derive d 
from observations of un-reddened quasars dHopkins et alj l2007t) . 
where the luminosity is determined using the accretion rate from 
the GADGET-2 simulations as described above. 

The dust distribution is determined by projecting the total gas- 
phase metal density in the GADGET-2 simulations on to a 3-D 
ada ptive mesh refinement grid ass uming a dust-to-metal ratio of 
0.4 dDwekll 19981 : 1 James et alj|2002h . We have used a maximum re- 
finement level of 10, which results in a minimum cell size of 55k" 1 
pc. This refinement level is sufficient to ensure the SEDs are con- 
verged to within a few per cent because the structure present in 
the GADGET-2 simulations is sufficiently resolved; if the resolu- 
tion of the GADGET-2 simulations were finer, the radiative trans- 
fer would require correspondingly smaller cell sizes. Note that we 
assume the ISM is smooth on s cales below the GA DGET- 2 resolu- 
tion and do not make use of the lGroves et al. liooi) sub-resolution 
photodissociation region model. The details of, motivation for, and 
implications of this choice are discussed in sections 2.2. 1 and 4.6 
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of HI 1 . We assume the d ust has properties given by the Milky Way 
R = 3.1 dust model of IWeingartner & Draind {2001 1) as updated 
bv lDraine & Li|j2007t) . The (sub)mm fluxes are similar if the LMC 
or SMC dust models are used. 

Once the star and BH particles are assigned SEDs and the dust 
density field is specified, SUNRISE performs the radiative transfer 
using a Monte Carlo approach by emitting photon packets that are 
scattered and absorbed by dust as they propagate through the ISM. 
The energy absorbed by dust is re-radiated in the IR. Dust tem- 
peratures, which depend on both grain size and the local radiation 
field, are calculated by assuming the dust is in thermal equilibrium. 
The ISM of our simulated galaxies can often be optically thick at 
IR wavelengths, so SUNRISE calculates the effects of dust self- 
absorption using an iterative method. This is crucial for ensuring 
accurate dust temperatures. 

The SUNRISE calculation yields spatially resolved SEDs 
(analogous to integral field unit spectrograph data) of the simulated 
galaxies viewed from different viewing angles. Here, we use 7 cam- 
eras distributed isotropically in solid angle. We use the SCUBA-2 
850- [im, AzTEC 1 . 1 -mm, and ALMA B ands 6 and 7 filter response 
curves to calculate the (sub)mm flux densities. Depending on the 
mass and the SMG flux cut, the simulations are selected as SMGs 
for ~ — 80 snapshots, and each is viewed from 7 viewing angles. 
Consequently, we have a sample of ~ 3.7 x 10 4 distinct synthetic 
SMG SEDs that we use to derive fitting functions for the isolated 
disc and galaxy-pair SMGs and duty cycles for the starburst SMGs. 



4 PREDICTING (SUB)MM NUMBER COUNTS 

To calculate the total SMG number counts predicted by our model, 
we must account for all subpopulations, including the infall-stage 
galaxy-pair SMGs discussed in Hll and H12, late-stage merger- 
induced starbursts, and isolated discs. To calculate the counts for 
the two subpopulations associated with mergers, we must combine 
the duty cycles of the mergers [the time the merger has (sub)mm 
flux greater than some flux cut] with merger rates because the 
number density is calculated by multiplying the duty cycles by the 
merger rates. For the isolated discs, we require the number density 
of a disc galaxy as a function of its properties and the (sub)mm flux 
associated with that galaxy. We describe our methods for predicting 
the counts of each subpopulation now. 

4.1 Late-stage merger-induced starbursts 

To predict the number counts for the population of late-stage 
merger-induced starburst SMGs, we combine merger rates - which 
depend on mass, mass ratio, gas fraction, and redshift - from the 
SEM with (sub)mm light curves from our simulations. We use the 
(sub)mm light curves from the simulations directly because it is dif- 
ficult to analytically model the dynamical evolution of the mergers, 
which can depend on the galaxy masses, merger mass ratio, progen- 
itor redshift, gas fraction, and orbital properties. For the SMG sub- 
population attributable to mergers, the number density of sources 
with flux density greater than S\ at redshift z is 



i(> S\, z) = 



dN{> S\,z) 



dV 



where dN/dVdtd\og Mb al dfidf g (M^ aY , /j,, f g , z) is the number 
of mergers per comoving volume element per unit time per dex 
baryonic mass per unit mass ratio per unit gas fraction, which is a 
function of progenitor baryonic mass A/bar, merger mass ratio (x, 
gas fraction at merger f g , and redshift z, and t(S\, Mbar, M> /s> z ) 
is the amount of time (duty cycle) for which a merger with most- 
massive-progenitor baryonic mass Aibar, mass ratio [i, and gas 
fraction f g at redshift z has flux density > S\. 

4.1.1 Duty cycles 

We calculate the duty cycles r(Ss5o) and r(Si.i) for various 5*850 
and Si.i values for the late-stage merger-induced starburst phase of 
our merger simulations. We neglect the dependence of duty cycle 
on gas fraction because sampling the range of initial gas fractions 
in addition to masses, mass ratios, and orbits is computationally 
prohibitive. Instead, as described above, we initialise the mergers 
with gas fraction f g — 0.8 so that sufficient gas remains at merger 
coalescence Q 

Similarly, because of computational limitations, we scale all 
initial disc galaxies to z ~ 3. We will see below that all else being 
equal, the dependence of (sub)mm flux density on z is small (< 
0.13 dex) for the redshift range of interest (z ~ 1 — 6), so we 
assume the duty cycles are independent of redshift and place the 
mergers at z = 3 (which is approximately the median redshift for 
our model SMGs) when calculating the duty cycles. Note, however, 
that the submm duty cycles for the starbursts may differ for mergers 
with progenitor disc properties scaled to different redshifts, but our 
model does not capture this effect. However, because most SMGs 
in our model have z ~ 2 — 4 (see below), this uncertainty should 
be subdominant. 

For each S\, we average the duty cycles for each set of mod- 
els with identical (-Mbar; A*) and men fit me resulting r (A/bar, 
surface with a second-degree polynomial in A/bar and /j, to estimate 
the duty cycle for (A/bar, y) values not explicitly sampled by our 
simulations. 



4.1.2 Merger rates 

The other ingredient needed to predict the counts for merger- 
induced starbursts is the merger rates. We use rates from the SEM 
described in detail in lHopkins et al. (2008c, 2010a. he), which we 
will briefly summarise here. The model starts with a halo mass 
function that has been calibrated using high-resolution TV-body 
simulations. Galaxies are assigned to haloes using an observed 
SMF for star-forming gala xies and the halo occupation formalism 
dConrov & W echsler 20091) . We use a fiducial SMF that is a com- 
bination of multiple observed SMFs, in which each covers a subset 
of the total redshift ra nge. For z < 2, w e use the SMF of star- 
forming galaxies fr omlllbert et alj d2010l). F or 2.0 ^ z ^ 3.75, 
we use the SMF of iMarchesini et al. I d2009l) because their survey 
is amongst the widest and deepest available and because they have 
performed a thorough analysis of the random and systematic un- 
cer tainties affecting the SMF determination. For z > 3.75, we use 
the lFontana et al.1 1 20061) SMF parameterization; though they only 



dN 



-(M hat ,fl,fg,z) 



dVdtd log A4bnrdfj,df g 
x T(S\,M baT ,ii,f g ,z)dlogM hsl , T d[idf g , (1) 



7 Note, however, that we calculate the starburst duty cycle using only the 
snapshots that sample the final starburst induced at merger coalescence, so 
the gas fraction is typically less than 40 per cent. We treat the early-stage 
galaxy-pair SMG contribution separately below. 
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Figure 1. Number density of disc galaxies, dN/ dV dlog(M* /M©) 
(Mpc~ 3 dex -1 ), versus M*(Mq) for integer redshifts in the range z = 
— 6 for our compo site SMF. For z < 2, we use the SMF for star- 
forming galaxies from l llbert et alj l2010l) . For 2 ^ z ( 3.75, we use the 
iMarchesini et alj J2009h SMF, and for z > 3.75, we use the Fontana et al. 
12006) parameterisation of the SMF. 



constrained the SMF out to z ~ 4, the extrapola tion agrees rea- 
sonabl y well with the 4 < z < 7 constraints from I Gonzalez et al] 
J201 lj), so this extrapolation is not unreasonable. Because the SMF 
at z > 4 is uncertain, it may be possible to constrain the SMF at 
those redshifts by using the observed SMG redshift distribution and 
relative contributions of the subpopulations; we discuss the s e pos- 
sibilities below. The interested reader should see iHavward Em%) 
for a detailed exploration of how the choice of SMF affects the 
predictions of our model. We do not correct for the passive galaxy 
fraction beyond 2 > 2, but this fraction is relati vely small at z ~ 2 
and decreases rapidly at higher redshifts (e.g.. IWuvts et alj|201 ll ; 
iBrammer et alj[201 lh . Our composite SMF at integer redshifts in 
the range z = — 6 is plotted in Fig. [Tj Finally, we use halo-halo 
merger rates from high-resolution iV-body simulations and trans- 
late to galaxy-galaxy merger rates by assuming the galaxies merge 
on a dynamical friction time-scale. 

The merger rates, SMF, and observed gas fractions are all un- 
certain. The merger rates are uncertain at the factor of ~ 2 level; the 
various sources of uncertainty and ef fects of modify i ng the model 
assumptions are discussed in detail in IHopkins et ail < l2010d lbl). At 
the redshifts of interest the random and systematic uncertainties in 
the SMF are comparable to the total uncertainty in the merger rates. 



4.1.3 Predicted counts 

Using the above assumptions, Equation (TJ becomes 

dN 



n(> S\, z) 



(Af bar , fi, z) 



dVdtd log Mbai-d^ 
x T(S\,M b a I ,[i)dlogM ba , I dv, 



(2) 



To calculate the observable cumulative counts (deg~ ), we 
must multiply by dV/dQdz, the comoving volume element in solid 
angle dQ and redshift interval dz, and integrate over redshift: 



dN(> Sx) _ f dN(> S x ,z) dV 



dn 
where 



dV dttdz 



(z)dz, 



(3) 



dV c (l + z) 2 D 2 A (z) 

-\ z ) — 



(4) 



dQdz y ~' H E(z) 
Here, Da (z) is the angular diameter distance at redshift z and 

4.2 Isolated discs and early-stage mergers 

We treat the isolated discs and early-stage mergers, which are dom- 
inated by quiescent star formation, in a semi-empirical manner, in 
which we assign galaxy properties based off observations. To cal- 
culate the observed (sub)mm flux densities using scaling relations 
similar to those of HI 1, we must determine the SFR and dust mass 
of a galaxy as a function of stellar mass and redshift. We then use 
SMF and merger rates to calculate the (sub)mm counts for these 
populations. 

4.2.1 Assigning galaxy properties 

Following IHopkins et al. I feOlOallcT) . we assign gas fractions and 
sizes as a function of stellar mass using observationally derived 
relations. We present the relevan t relations below, but we refer the 
reader to IHopkins et al. for full details, including the 

list of observati ons used to derive the r elations and justifications for 
the forms used. IHopkins ai] d2010d) have shown that this model 
reproduces global constraints, such as the IR luminosity function 
at various redshifts and the star formation history of the Universe, 
among others; these results support the application of the model in 
this work. 

The baryonic gas fraction, / gas = M gas /(M gaa + M*), of a 
galaxy of stellar mas s M* and redshift z, a s determined from the 
obser vations listed in IHopkins et al" N2010ch , is given by Equation 
(1) of lHopkins etatl feOlOch ~ 



/ gas (M*|z = 0) 
/ gaB (M*,z) 



fo 



= fo 



1 + (M*/10 915 M ) - 4 ' 

l-T(*)[l-/ S/aNV 



(5) 



where r(z) is the fractional look-back time to redshift z. At a given 
mass, galaxy gas fractions increase with redshift. At fixed redshift, 
they decrease with stellar mass. Using / gas (M*, z), we can calcu- 
late the gas mass as a function of Ai* and 2, 



M gas (M t ,2) 



/ gaB (M*,2) 
l-/ga S (M»,2) 



(6) 



Similarly, we parameterize the radius of the gas disc as a 
function of mass and redshift using the observations listed in 
IHopkins et al. (2010c). (Note that the stella r disc radii are signif- 
icantly smaller.) The relation ( Equation 2 of lHopkins et alj|2010d ; 
see also lSomerville et"al]|2008h is 



R e {lvL\z = 0) = Ro = 5.28 kpc 
R^M^z) = i? (l + 2)" 0,6 



AL 



10 10 M s 



(7) 



(8) 



We assume the quiescently star-fo rming discs obey the KS re- 



We assume the quiescently star-to 
lation jSchmidtll 1 9591 : iKennicutdl 1 998h 



= 1.3 x 10" 4 M yr -1 kpc- 



M oP c- 



(9) 



where E* and E gas are the SFR and gas sur face densities, re- 
spectively and 71k = 1-4 dKennicuttl 1 19981) , at all redshifts, 
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Figure 2. SFR ( M yr' 1 ) versus M*( Mq) for model disc galaxies at 
integer redshifts in the range z = — 6 (solid lines) and from the observa- 
tionally derived fitting function of Whitaker et al] J20 12b at z = 0, 1, and 
2 (dashed lines). The normalisation of the relation increases with redshift 
both because gas fractions are higher and galaxies are more compact. The 
model agrees reasonably well with the observations except at z ~ 0, but 
we shall see that the z ~ contribution to ou r model SMG population is 
small. If the SFRs from I Whitaker et al J <2012l) were used instead of those 
calculated from Equation ilO\ . the z ~ contribution would be even less, 
so the discrepancy is unimportant. 



as is supported by observations (e . g.. [Narayanan et alj |2008| : 



iDaddi et alj|201fj: iGenzel et alj[201ol ; iNaravanan et al 



,f. iNaravanan et alj l2012bf) T We normalise the relation 



2011 



but 



assum- 
ing a iKroupal j200ll) IMF. Assuming S gas ss M gaii /(nRl) and 



M*(M*,z) = 



where M* is the SFR, we find 

10 4 V K_1 ( M Ha8 (M*,«) 



1.3 



kpc 



-2(r, 



10 10 Mq 
-1) 

M yr 



(10) 



which can be recast in terms of M* rather than M gas using Equa- 
tions (|5j and J6). Fig.|2]shows the SFR- Ah relation given by Equa- 
tion l |10b for integer r edshifts in the range z = — 6 and the ob- 
served relations from IWhitaker et al.1 j2012l) for z ~ 0, 1 , and 2. 
The agreement between our relation and those observed for z ~ 1 
and 2 is reasonable for the masses (M* > 10 11 M Q ) relevant to 
the SMG population. The agreement is less good for z ~ 0, but, 
as we shall see below, this is u nimportant because th e fraction of 
SMGs at z < 1 is small. If the I Whitaker et alj d2012l) SFRs were 
used instead of those calculated from Equation JlOt . the z ~ 
contribution would be even smaller. 

In addition to the SFR, we require the dust mass to calcu- 
late the (sub)mm flux densities. To determine the dust mass, we 
must know the gas-phase metallicity. Observations have demon- 
strated that metallicity increases with stellar mass; t his relationship 



has been constrained for redshifts z 







3.5 1 Tremonti et al 
1 1 — 



20041: Savaglio et al. 2005; lErb et al. 2006; Kewl ev & Ellisor 



20081 : iMaiolino et alj|2008h . iMaiolino et all J2008I) parameterized 
the evolution of the mass-metallicity relation (MMR) with redshift 
using the form 

12 + log(0/H) = -0.0864 [log A/* -logM (z)] 2 

+ A'o(z). (11) 
They determine the values of log Mo and Ko at redshifts z = 



0.07, . 7, 2.2, and 3.5 using th e observations of l Kewlev & Ellison! 
j2008h . ISavaglio et all d2005l) . lErb et al.1 J2006h . and their own 
work, respectively. To crudely capture the evolution of the MMR 
with redshift, we fi t the v alues of log Mo and Ko given in Table 
5 of IMaiolino et al.1 d2008h as power laws in (1 + z); the result is 
logMo(z) « 11.07(1 + z) 094 and K (z) 9.09(1 + 



Using 12 + log(0/H) G = 8.69 dAsplund et a"i]|2009l) . we 



have 

log(0/H) 



- log(O/H) = 

- 0.0864 [log Af* - 11.07(1 + z) 
+ 9.09(1 + z)" 017 - 8.69. 



0.94] 2 



(12) 

The solar metal fraction is Z Q = 0.0142 ( lAsplund et alj2009l) . so 
Z(M„, z) = 0.0142 (io 1 °g(°/ H )- 1 °g( o / H )0) . (13) 

We assume the dust mass is proportional to the gas-phase metal 
. Thus, 



M d {M t ,z) = hU 



- fga S (M*,z) 



X Z(M„,z)f dt n 



(14) 



where we use a dust-to-metal ratio /dtm = 0.4 l lDwe' 
Ijames et alj|2002h . 

Motivated by equation (1) of Hll, we fit the (sub)mm flux 
densities of our simulated galaxies as power laws in SFR and A/d. 
We find that 



S850 



0.81 mJy 

M* 



100 M yr- 1 
S1.1 = 0.35 mJy 

100 M Q yr" 1 



Af d 



10 s M 



10 8 M 



(15) 



(16) 



is accurate to within 0.13 dex for z ~ 1 — 6. (The flux for galaxies 
at z < 0.5 is underestimated significantly by these equations, but 
such galaxies contribute little to the overall counts because of the 
smaller cosmological volume probed and the significantly lower 
gas fractions and SFRs, so this underestimate is unimportant for 
our results.) The (sub)mm flux is insensitive to redshift in this red- 
shift range because as redshift increases, the decrease in flux caused 
by the increased luminosity distance is almost exactly cancelled by 
the increase in flux caused by the rest-frame wavelength moving 
closer to the peak of the dust emiss ion (this effect is referred to as 
the negative K correction; see, e.g. jBlain et "aT] |2002l) . By combin- 
ing Equations l |10t and J 1 4t with Equations \\5\ and (QJo), we can 
calculate Ss5o(M*, z) and 5i.i(M*, z), 



5'sr, 



0.81 mJy 



0.013 



AL 



10 4 \°V M e; 



10 10 M 



Re 

kpc 



Si. 



10 8 M 

0.35 mJy 



0.013 



M s 



10 4 



Z(M*,z)/dt„ 



( M g; 



(17) 



^lO 10 Mq 



Re 

kpc 



10 s M 



Z(M*, z)/dt„ 



(18) 
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Figure 3. Observed-frame 850-|j.m (5g5o; top) and 1.1-mm (Si.i; bottom) 
flux density in mJy versus M*( Mq) for isolated discs at integer redshifts 
in the range z = 1 — 6 (see Equations 1 171 and 118) . The (sub)mm flux of 
a disc of fixed iW* increases with redshift for two reasons: 1. As shown in 
Fig. [5] the normalisation of the SFR-Af* relation increases with redshift. 2. 
For fixed M* , gas fraction increases with redshift. For fixed Z, a higher gas 
fraction corresponds to a higher gas-phase metal mass. However, because 
the normalisation of the MMR decreases as z increases, the increase of 
the gas-phase metal mass with gas fraction is partially mitigated. Both the 
increased SFR and increased dust mass cause the (sub)mm flux to increase. 

where we can substitute the appropriate expressions for M gas , R e , 
and Z to express Ssso and Si.i in terms of M* and z only. 

Fig.[3]shows the Ssso — M* and S\,\ — M* relations given by 
Equations (TTJl and ( 1 1 3b , respectively, for isolated discs at integer 
redshifts in the range z — 1 — 6. As redshift increases, galaxies 
become more gas-rich and compact; both effects cause the SFR 
for a given A/* to increase (see Fig. |2](. For fixed Z, a higher gas 
fraction corresponds to a higher gas-phase metal mass. However, 
because the normalisation of the MMR decreases as z increases, 
the increase of the gas-phase metal mass with gas fraction is par- 
tially mitigated[f|The increased SFR and Md both result in a higher 
(sub)mm flux for a given M*. To produce an isolated disc SMG 
(•Ssso ^3 — 5 mJy, or Si.i > 1 — 2 mJy) at z ~ 2 — 3, we re- 

8 For very high redshifts, at which there should be very few metals and thus 
very little dust, the (sub)mm flux of a galaxy of a given mass should de- 
crease sharply. However, the precise redshift dependence depends on still- 
uncertain details of dust production; thus, it is possible that constraints on 
the fraction of SMGs with z > > 4 will yield insight into the physics of 
dust production. 



quire M+ > 10 11 Mq. This va lue is consistent with the results of 
iMichalowski et alj J20ldl2012h . 

Note also that we can use these relations to calculate the ex- 
pected Ssho/Sx.i ratio , Ssso/ Si. i ~ 2.3 (this is similar to obser- 
vational estimates; e.g. jAustermann et alj201(j|) . For simplicity, we 
will use this ratio to derive approximate Ssso values to also show 
an Ssso axis on the relevant plots. 

4.2.2 Inf all-stage galaxy-pair SMGs 

During the infall stage of a merger, the discs are dominated by 
quiescent star formation that would occur even if they were not 
merging. Only for nuclear separation < 10 kpijf] do the discs have 
SFRs that are significantly elevated by the mutual tidal interactions 
(HI 2); this re sult is consistent w ith observed SFR elevations in 
mergers (e.g., IScudder et alJl2012T) . Thus, during the infall stage, 
we assume the discs are in a steady state (i.e., they have con- 
stant SFR and dust mass); even without a source of additional gas, 
this is a reasonable approximation for the infall stage to within a 
factor of < 2 (see fig. 1 of Hll). For a merger of two progeni- 
tors with stellar masses M*,i and M*,2, the total flux density is 
S\ = SxiM+,1) + Sx(M*,2). The typical beam sizes of single- 
dish (sub)mm telescopes are 15", or ~ 130 kpc at z ~ 2 — 3; 
schematically, when the projected separation is less than this dis- 
tance, the sources are blended into a single source0 To predict 
single-dish counts, we assume the galaxies should be treated as a 
single source if the physical separation is < 100 kpc. From our 
simulations, which use cosmologically motivated orbits, we find 
that this time-scale is of order ~ 500 Myr. Though the time-scale 
depends slightly on the most-massive-progenitor mass, we neglect 
this dependence because it is subdominant to various other uncer- 
tainties. However, this time-scale is derived from the z ~ 2 — 3 
simulations and thus may be too long for mergers at higher z. Given 
the above assumptions, the duty cycle for a given S' x and merger 
described by more-massive progenitor mass M*,i and stellar mass 
ration = M*, a /M*,i is 0.5 Gyrif Sa(M*,i) + 5x(M*, 1 /i) > S' x 
and otherwise. With the duty cycle in hand, we can use Equations 
([2J and ([3j to calculate the predicted number density and counts. 

To predict counts for ALMA, we simply assume that the two 
discs are resolved into individual sources and thus treat them as two 
isolated disc galaxies, as described below. 

4.2.3 Isolated disc counts 

For a given S\ and z, we invert the 5a (M* , z) functions (Equations 
[T71 and |ISt to calculate the minimum M* required for a galaxy 
at redshift z to have (sub)mm flux density > S\, M+(S\\z). To 
calculate the number density n(> S\, z), we then simply use the 
star-forming galaxy SMF to calculate 

n{> Sx,z) =n[> AL{Sx\z),z], (19) 

9 This value is derived for the z ~ 3 simulations presented here, and it may 
differ for z ~ simulations because of differences in structural properties 
and gas fractions. 

10 Because of the computational expense involved, we do not create syn- 
thetic maps, add noise, convolve the maps with a Gaussian beam, and then 
calculate the fluxes of the sources. However, the uncertainty caused by our 
simple method of calculating the total flux is subdominant to other uncer- 
tainties inherent in the model. 
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Figure 4. Predicted cumulative number counts for the unlensed SMG pop- 
ulation as observed with single-dish (sub)mm telescopes, N(> Si.i), 
in deg~ 2 , versus Si.i (mJy). The counts are decomposed into the three 
SMG subpopulations we model: the green long-dashed line corresponds to 
isolated disc galaxies, the blue dashed to galaxy-pair SMGs (i.e., infall- 
stage pre-starburst mergers), and the red dash-dotted to merger-induced 
starbursts. The black solid line is the total for all SMG subpopulations 
we model. The mo del predictions of iBaugh et alj fc005j, navy dotted), 
iGranato et alj j2004l magenta dotted), and iFontanot et al J 12007L maroon 
dotted) are shown for comparison. The points are observed 1.1 -mm and 
850- and 870-p.m counts (see the text for details). The 850- and 870-^rn 
counts have been converted to 1.1 -mm counts by assuming S&50/S1.1 ~ 
S'87o/>Si.i = 2.3. This ratio has also been used to show the approxi- 
mate Sgso on the top axis. The hatched area shows the regime where weak 
lensing is expected t o significantly boost the counts for overdense fields 
lAretxaga et al.|[2oTlh . The counts predicted by our model agree very well 
with the counts that are not thought to be boosted significantly by lensing. 
N.B. The steepness of the cutoff in the starburst counts at S1.1 > 4 mJy is 
artificial; see the text for details. 



and we use Equation ® to calculate the predicted counts. To pre- 
dict counts for single-dish (sub)mm telescopes, where the galaxy- 
pair SMGs are blended into a single source, we subtract the fraction 
of galaxies with M* > M ir (Sx\z) that are in mergers from the iso- 
lated disc counts to avoid double counting. 



5 RESULTS 

Here, we present the key results of this work, the SMG cumulative 
number counts, the relative contributions of the subpopulations, 
and the reds hift distribution pr edicted by our model. We focus on 
the AzTEC (Wi lson et alj200ct) 1.1-mm counts here because to our 
knowledge, the best-constrained blank-field counts (i.e., those from 
the deepest and widest surveys) have been determined us ing that 
instrument ( lAustermann et al]|201fj| : lAretxaga et alj|201lh . How- 
ever, because for our simulated SMGs, Sg$o/Si.i = 2.3 to within 
~ 30 per cent, we can easily convert the 1.1-mm counts to 850-|j.m 
counts. Thus, we include both S1.1 and Ssso values on the relevant 
plots and convert observed 850- vim counts to 1.1-mm counts by 
assuming the same ratio holds for real SMGs. 



5.1 SMG number counts for single-dish observations 

Fig. E] shows the total cumulative 1.1-mm number counts (black 
solid line), which are calculated from the cumulative number den- 
sity using Equation {3). We decompose the counts into isolated 



discs (green long-dashed), galaxy pairs (blue dashed), and star- 
bursts induced at merger coalescence (red dash-dotted); the relative 
contribution of each subpopulation is discussed in Section [5~2l The 
data points in Fig.[4] are observed count s from va rious surveys: 1.1 



mm counts from 



Aretxagaetal.l l 2011 



201(1 squares). lHatsukade et all d201 1 



circles), lAustermann et al 



diamonds), and Scott et al 



20 id . trian gles); 850-yim count s from iKnudsen et alj I2OO8L as 
teris ks) andlZemcov et alj J2010t plus signs); and 870- (im counts 
from IWeifi et alj d2009l . x's). The 850- and 870-iim counts have 
been converted to 1.1-mm counts by assu ming Sgso/Si.i ~ 
SWo/ffi.i = 2.3. T he mo del predictions of Baugh et"al] J2005h . 
IGranato et alj d2004l) . and IFontanot et alj d2007t> are shown for 
comparison. 

The predicted and observed counts are in good agreement at 
the lowest fluxes, but the predicted counts are less than som e of 
those observed at the bright end. The lAustermann et alj j201Ch and 
lAretxaga et alj J201 lh surveys are the two largest (~ 0.7 deg -2 ), 
so their counts should be least affected by cosmic variance and thus 
most robust. Thus, it is encoura ging that the agreement b etween 
our predicted counts and those of lAustermann et all l l2010l) is very 
good at all fluxes. Th e disagreement betwee n our predicted counts 
and those observed bv lAretxaga et alj j201 ll) is significant even for 
the lower flux bins (a factor o f ~ 2 for the S1.1 > 2 mJy bin). 
However. lAretxaga et al.11201 ll) conclude that the excess of sources 
at S1.1 > 5 mJy compa red with the SHADES field observed by 
lAustermann et alj J2010h is caused by sources moderately ampli- 
fied by galaxy-galaxy and galaxy-group lens ing. At higher fluxes , 
the effect of lensing is more sign ificant jNegrello et alj 1 20071 : 
IPaciga etalj|2009l ; iLima et alj[201CT) . and it would be incredibly 
difficult to e xplain the sources with mm flux density >> 10 mJy 
observed by I Vieira et ail J2010h and lNegrello et alj d2010h if they 
are not strongly lensed. We do not include the effects of gravita- 
tional lensing in our model, so it i s unsurprising t hat w e signif- 
icantly under- predict the counts of lAretxaga et alj J201ll) despite 
the excellent agreement between our counts and those observed by 
lAustermann et alj d20 10h . 

Additionally, it is important to note that the steepness of the 
cutoff in the starburst counts at S1.1 > 4 mJy is artificial: because 
we determine the fluxes of the isolated discs and galaxy pairs in 
an analytic way, we can extrapolate to arbitrarily high masses for 
those populations. For the starbursts however, we are limited by the 
parameter space spanned by our merger simulations. None of our 
starburst SMGs reach S1.1 ^ 6.5 mJy (or Ssso <; 15 mJy), so the 
duty cycle for all starbursts for S1.1 6.5 mJy is zero. If we were 
to simulate a galaxy more massive than our most massive model 
(b6), the simulation would reach a correspondingly higher flux, so 
the predicted counts for S1.1 6.5 mJy would no longer be zero. 
However, the rarity of such objects does not justify the additional 
computational expense. Thus for S1.1 > 4 mJy (or Ssso > 9 mJy), 
the starburst counts should be considered a lower limit. A simple 
extrapolation from the lower-flux starburst counts suggests that our 
model may even over-predict the counts of the brightest sources. 
However, the observed number density of sources with S1.1 > 5 
mJy is highly uncertain because of the effects of small number 
statistics, cosmic variance, and lensing, and the uncertainty in the 
model prediction is significant because of uncertainties in the abun- 
dances and merger rates of such extreme systems. Thus, the counts 
for the brightest sources should be interpreted with caution. 

Furthermore, we do not attempt to model some other poten- 
tial contributions to the SMG population. In particular, we do not 
include contributions from mergers of more than two discs, clus- 
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Table 4. Single-dish-detected SMG cumulative number counts 



11 



S-l.x" ~ Ssbo b N(> Si.iY Fractional contribution 



(mJy) 


(mJy) 


(deg- 2 ) 


Isolated discs'' 


Galaxy pairs c 


Starbursts' 


U.J 


i i 
t . t 




V.JZ 


U. jj 


r\ 1 t 
U. 1 j 


1 .U 


L.j 


ZD 10 


n 

U.zj 


U.36 


U.39 


1 .J 


5. 4 


nil 


U.2U 




U.Zo 


2.0 


4.6 


354 


0.14 


0.51 


0.35 


2.5 


5.7 


200 


0.09 


0.45 


0.46 


3.0 


6.9 


110 


0.06 


0.41 


0.52 


3.5 


8.0 


65 


0.04 


0.38 


0.58 


4.0 


9.2 


41 


0.03 


0.33 


0.64 


4.5 


10.3 


27 


0.02 


0.28 


0.70 


5.0 


11.5 


13 


0.01 


0.29 


0.69 


5.5 


12.6 


6 


0.01 


0.40 


0.59 


6.0 


13.8 


2 


0.01 


0.53 


0.45 



" 1.1-mm flux density. * Approximate 850-i^m flux density calculated assuming Ssso/Si.i = 
2.3. c Cumulative number counts of SMGs with 1.1-mm flux density greater than the S\.\ value 
given in the first column. ' Fractional contribution of each subpopulation to the total cumulative 
counts for the given Si.i. 



2.3S„ = S 850 (mJy) 



10 



1.0 



4 0.8 



o 0.6 



0.4 



■g 0.2 



0.0 



isolated discs 
galaxy pairs 
starbursts 



S„ (mJy) 

Figure 5. Fractional contribution of each subpopulation to the total cumu- 
lative counts versus Si.i. The lines are the same as in Fig. [4] At the lowest 
fluxes, the isolated discs dominate, whereas at the highest fluxes, the star- 
bursts dominate. The galaxy pairs are ~ 30-50 per cent of the population at 
all fluxes plotted here. 



ters, or phy sically unrelated s ources blended into a single (sub)mm 
source (see lWang et al.|[201 ll for evidence of the last type). 

Given these caveats and the modelling uncertainties, our pre- 
dicted counts are clearly consistent with those observed, and in- 
cluding lensing and the previously mentioned additional possible 
contributions to the SMG population would tend to increase the 
number counts. Also, we stress that our model is conservative in 
the sense that it uses a Kroupa - rather than top-heavy or flat - IMF 
and is tied to observations whenever possible. The consistency of 
the predicted and observed counts suggests that the observed SMG 
counts may not provide evidence for IMF variation; this will be 
discussed in detail in Section[6~T| 



5.2 Relative contributions of the subpopulations 

In previous work (H 1 1 ; lHavwardetalj2011bl ; H 12), we argued that 
the SMG population is not exclusively late-stage merger-induced 
starbursts but rather a heterogeneous collection of starbursts, infall- 



stage mergers ('galaxy-pair SMGs'), and isolated discs. However, 
so far we have only presented the physical reasons one should ex- 
pect such heterogeneity. It is crucial to quantify the relative impor- 
tance of each subpopulation, so we do this now. 

The counts shown in Fig. [4] are divided into subpopulations, 
but the relative contributions can be read more easily from Fig. [5] 
which shows the fractional contribution of each subpopulation to 
the total cumulative counts. At the lowest fluxes, the isolated disc 
contribution is the most significant. At Si.i ~ 0.8 mJy (Ssso ~ 2 
mJy), the three subpopulations contribute almost equally. As ex- 
pected from conventional wisdom, the starbursts dominate at the 
highest fluxes. However, contrary to conventional wisdom, the 
bright SMGs are not exclusively merger-induced starbursts: from 
Fig. [5] we see that at all fluxes plotted, the galaxy pairs account for 
~ 30-50 per cent of the total predicted counts, so they are a signif- 
icant subpopulation of SMGs in our model. As explained in HI 2, 
the galaxy-pair SMGs are not physically analogous to the merger- 
induced starburst SMGs; thus, their potentially significant contribu- 
tion to the SMG population can complicate physical interpretation 
of the observed properties of SMGs. 

It is interesting to compare the relative contributions of the 
isolated disc and galaxy-pair subpopulations because the relative 
contributions can be understood - at least schematically - in a sim- 
ple manner. For a major merger of two galaxies with = Mi so , 
the flux of the resulting galaxy-pair SMG is approximately twice 
that of the individual isolated discs, 2Si.i (Mi so ). Because Si.i de- 
pends sub-linearly on (see Fig. [3}, for an isolated disc to have 
jSi.1 equal to that of the galaxy pair, it must have > 3Mi so . 
Thus, the relative contribution of the two subpopulations depends 
on whether the number density of M* = 3Afi so discs divided by 
that of A4* = Mi so discs, n(3Afi so )/n(Afi so ), is greater than the 
fraction of Af* = Afi so discs undergoing a major merger, which 
is the merger rate times the duty cycle of the infall phase (~500 
Myr). If the former is larger, the Af* = 3Afi so discs will dominate 
the pairs of M+ = Afi ao discs, whereas if the merger fraction is 
higher than the relative number density, the galaxy pairs will dom- 
inate. 

The latter scenario is likely for bright SMGs, which are on 
the exponential tail of the SMF. For example, at z ~ 2 — 3, 
a galaxy with Af* = 10 11 Mq undergoes ~ 0.3 mergers per 
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Gyr. Thus, if we assume a duty cycle of 500 Myr for the galaxy- 
pair phase, approxim ately 15 per cent of such galaxies will be in 
galaxy pairs. For the iMarchesini et alj d2009h SMF, the number 
density of Af* — 3 x 10 11 Mq galaxies is ~ 8 per cent that of 



10 1 



Mq galaxies. Therefore, by the above logic, the pairs 
in 



Af* = 

of Af* = 10 11 Mq galaxies will contribute more to the submm 
counts than the isolated M* = 3 x 10 11 Mq discs. This simple ar- 
gument demonstrates why the galaxy pairs become dominant over 



the isolated discs for Si. i > 0.7 mJy (Sg 



> 



1.6 mJy). However, 



the threshold for dominance depends on both the Si.i — Ah scaling 
and the shape of the SMF at the high-mass end. Thus, observation- 
ally constraining the fraction of the SMG population that is galaxy 
pairs can provide useful constraints on both the (sub)mm flux-Af* 
relation and the shape of the massive end of the SMF. 

Unfortunately, the relative contribution of the starburst sub- 
population cannot be explained in as simple a manner. The duty 
cycles for the merger-induced starbursts depend sensitively on pro- 
genitor mass and merger mass ratio, so the mapping from merger 
rate to number density is not as simple as it is for the isolated discs 
and galaxy pairs. Fortunately, the SMF uncertainty, which is very 
significant for the overall counts, is relatively unimportant for the 
relative contribution of starbursts and galaxy pairs. Thus, the rel- 
ative contributions of starbursts and galaxy pairs depend primar- 
ily on their relative duty cycles. (To achieve a given flux density, 
one requires a less massive starburst than galaxy pair because the 
starburst increases the (sub)mm flux density moderately. Thus, the 
relative number density also matters. However, the inefficiency of 
starbursts at increasing the (sub)mm flux density of the system pre- 
vents significantly less massive but more common starbursts from 
dominating over more massive and rarer galaxy pairs.) The duty 
cycles are uncertain, but given that in our fiducial model the galaxy 
pairs contribute ~ 30 — 50 per cent of the total counts and the un- 
certainty in the duty cycles is definitely less than a factor of 2 — 3, 
the prediction that both the starburst and galaxy pair subpopulations 
are significant (i.e., more than a few per cent of the population) is 
robust. 

Though there have been many observational hints suggest- 
ing the importance of the galaxy-pair contribution (see Hll and 
H12 for discussion), the physical importance of this subpopula- 
tion has to date not been fully appreciated, and the fractional con- 
tribution of galaxy-pair SMGs to the total counts remains rela- 
tively poorly constrained. However, clear observational evidence 
supporting the significance of th is subpopulation i s accumulat- 
ing: of the 12 SMGs presented in lEngel et al] bold) . 5 have CO 
emission that is resolved into two components with kinematics 
consistent with two merging discs. In two of the cases, the pro- 
jected separation of the two components is > 20 kpc; such ob- 
jects are prime examples of the galaxy-pair subpopulation. (See 
also lTacconi et all200dl2008tlBothwell et alj201CtlRiechers et alj 



1201 lal lbh lSmolcic et alj feOra) presented a larger sample of SMGs 
with (sub)mm-interferometric detections. They found that when 
observed with interferometers with < 2" resolution, ~ 15 — 40 
per cent of single-dish SMGs were resolved into multiple sources, 
which is consistent with our prediction for the relative contribution 
of the galaxy-pair subpopulation. ALMA observations will signif- 
icantly increase the number of SMGs observed with ~ 0.5" reso- 
lution and thus better constrain the galaxy-pair contribution to the 
SMG population. 

Further evidence for a galaxy-pair contribution consistent with 
what we predict is the fraction of the SMGs with multiple counter- 
parts at other wavelengths. One of the earliest observational indi- 
cations of this population came from the SCUBA 8-mJy survey: of 
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Figure 6. Top: the predicted redshift distribution of 1.1 -mm sources with 
Si.i > 1.5 mJy (Sg5o > 3.5 mJy; black solid line) compared with the 
observed distribution from lYun et alj j2012h (blue dashed line). The mean 
redshifts are 3.0 and 2.6 for the model and observed SMGs, respectively. 
Bottom: the predicted redshift distribution for sources with Si.i > 4 mJy 
(S85 5*, 9 mJy; black sol id line) compared with the observed distribution 
from ISmolcic et all j2012l) (blue dashed line). The mean redshifts are 3.5 
and 3.1 for the model and observed SMGs, respectively. In both panels, the 
flux limits were chosen to approximately match those of the observations. 
The brighter SMGs tend to be at higher redshifts. 



this sample of 850- ^tm sources. Ilvison et aLI J2002I) found that ~ 25 
per cent have multiple radio co unterparts. Appro ximately ten per 
cent of the GOODS-N 8 50-nm jPope et alj2006h. GOODS-N 1.1 



mm dChapin et alj|2009l) . SHADES 850- urn jlyison et al 
IClements etal.ll2008h . and GOODS-S 1.1 -mm jYun et al. 



20Q-/1 ; 
2012b 



sources have multiple counterparts. These fractions are somewhat 
smaller than the ~ 30-50 per cent contribution shown in Fig. 
f5] but both the predicted and observed fractions are uncertain. 
As explained above, the predicted fraction depends sensitively on 
the shape of the upper-end of the SMF and the relation between 
(sub)mm flux and Af*. Observations, on the other hand, may miss 
the more widely separated counterparts and cases when one of the 
counterparts is significantly more obscured or is radio-quiet. 



5.3 Redshift distribution 

In addition to the number counts, a successful model for the SMG 
population must reproduce the redshift distribution. Fig. [6] shows 
the redshift distribution of 1.1-mm sources predicted by our model 
for different 1.1-mm flux cuts (S1.1 > 1.5 mly, or Sgso > 3.5 
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Table 5. Single-dish-detected SMG redshift 
distribution 
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mJy, in the top panel and Si.i > 4 mJy, or Ssso > 9 mJy, in 
the bottom) along with some observed distributions that have sim- 
ilar flux limits. The redshift distributions are relatively broad, and 
they peak in the range z ~ 2 — 4 and decline at lower and higher 
redshifts. The Si.i > 1.5 mJy sources have mean redshift 3.0, 
whereas the Si.i > 4 mJy sources have mean redshift 3.5, so there 
is a tendency for the brighter s ources to be at higher redshifts; this 
trend agrees with ob servations dlvison et all2002t lYun et alj|2012l ; 
ISmolcic et al]|2012h . 

For the Si.i > 1.5 mJy SMGs (top panel of Fig.|6j, compared 
with the observations, our model predicts a higher mean redshift 
and a greater fraction of SMGs at z ~ 3 — 4. This discrep ancy may 
suggest that the extrapolation of the lFontana et al.l d2006h SMF we 
use for z > 3.75 over-predicts the number of massive galaxies 
at the highest redshifts. Furthermore, merger time-scales may be 
shorter at high redshift, and the dust content may be lower than in 
our models; both of these effects would decrease the high-redshift 
contribution. Additionally, constraining the redshift distribution is 
complicated by selection effects, counterpart identification, and, in 
some cases, a lack of spectroscopic redshifts. Typically, the selec- 
tion effects and counterpart identification make identifying higher- 
redshift SMGs more difficult. Fin ally, the significant differences 
amongst observed distributions (see lSmolcic et alj2012l for discus- 
sion) demonstrate the difficulty of determining the redshift distri- 
bution. Thus, the differences in the distributions may not be sig- 
nificant; spectroscopic follow-up of ALMA sources should clarify 
this issue (but note that the redshift distribution of ALMA sources, 
which is shown in the next section, should differ slightly). The typ- 
ical redshift for the brighter SMGs (bottom panel of Fig. [6} is also 
somewhat higher than that observed (3.5 compared with 3.1), and 
the peak at z ~ 1 — 2 is not reproduced. However, the observed 
distribution is based on a sample of 17 sources, so small-number 
statistics might explain the differences. 



5.4 Predicted ALMA-detected SMG number counts and 
redshift distribution 

In the above sections, we have presented our model predictions for 
surveys performed with single-dish (sub)mm telescopes such as the 
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Figure 7. Predicted cumulative number counts for the SMG subpopulation 
observed with telescopes, such as ALMA, that have resolution sufficient 
to resolve the galaxy-pair subpopulation into multiple sources. The green 
long-dashed line corresponds to isolated disc galaxies, the dash-dotted red 
to starbursts, and the black solid to the total counts. The single-dish counts 
from Fig.|4]are shown for comparison (blue dash dot dot dot). The top axis 
shows approximate Sgso values. Because some of the bright sources would 
be resolved into multiple sources with ALMA, the counts of bright ALMA 
sources are as much as a factor of 2 less than the single-dish counts. 
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Figure 8. Predicted redshift distributions for ALMA-detected SMGs with 
Si.i > 1.5 mJy (5850 ^ 3.5 mJy; black solid line; z = 3.0) and Si.i > 4 
mJy (Sgso > 9 mJy; blue dashed line; z = 3.5). As for the single-dish 
sources, the brighter SMGs tend to lie at higher redshift. 



JCMT, for which the ~ 15 arcsec beam causes the galaxy-pair 
SMGs to be blended into a single source for much of their evo- 
lution. If larger single-dish (sub)mm telescopes such as the LMT 
and CCAT are used, then the galaxy-pair SMGs will be blended 
for a smaller fraction of the infall stage of the merger. Thus, the 
duty cycle for the galaxy-pair phase will be shorter and the number 
counts of the bright sources consequently less. For interferometers 
with angular resolution < 0.5 arcsec (or ~ 4 kpc at z ~ 2 — 3), 
such as ALMA, essentially all galaxy-pair SMGs will be resolved 
into multiple sources because for such separations, the mergers are 
typically dominated by the starburst mode (H12). Thus, the galaxy 
pairs would contribute to the number counts as two quiescently 
star-forming galaxies rather than one blended source. 

To predict the counts and redshift distribution of ALMA 
sources, we modify our model by removing the galaxy-pair con- 
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Table 6. ALMA-detected SMG cumulative number counts 





~ ^850 


~ Salma-6' 


~ SaLMA-7 < ' 




Fractional contribution 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(deg~ 2 ) 


Isolated discs 


Starbursts^ 


U.J 


■ i 
i . i 


U.4 


n q 
U. o 


DOV / 


0.85 


0.15 


1 .U 


A.D 


U.o 


1 A 
1 .0 


ZU/Z 


0.50 


0.50 


1 .J 


J. 4 




-1 
Z.4 


4oZ 


0.56 


0.44 


z.u 


4.0 


1 .0 


J.Z 


1 A 
Z14 


0.41 


0.59 


2.5 


5.7 


2.0 


4.0 


126 


0.27 


0.73 


3.0 


6.9 


2.4 


4.8 


71 


0.19 


0.81 


3.5 


8.0 


2.8 


5.6 


43 


0.14 


0.86 


4.0 


9.2 


3.2 


6.4 


29 


0.09 


0.91 


4.5 


10.3 


3.6 


7.2 


20 


0.06 


0.94 


5.0 


11.5 


4.0 


8.0 


10 


0.05 


0.95 


5.5 


12.6 


4.4 


8.8 


3 


0.06 


0.94 


6.0 


13.8 


4.8 


9.6 


1 


0.08 


0.92 


" 1.1-mm 


flux density. 


b d Approximate 850- 


lim and ALMA Bands 6 & 7 flux densities calculated 


using the conver- 



sion factors Ssso/Si.i = 2.3, Salman/Si. l = and SalmA-7/>Si.i = 1-6. e Cumulative number counts 
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Table 7. ALMA-detected SMG redshift dis- 
tribution 



a 

Z 


iVinbi n /iVtotal(> Sl.l) 

1 mJy* 4 mJy'' 


0.0-0.5 


0.002 


0.006 


0.5-1.0 


0.021 


0.031 


1.0-1.5 


0.040 


0.047 


1.5-2.0 


0.116 


0.092 


2.0-2.5 


0.161 


0.094 


2.5-3.0 


0.257 


0.125 


3.0-3.5 


0.286 


0.224 


3.5-4.0 


0.043 


0.102 


4.0-4.5 


0.031 


0.094 


4.5-5.0 


0.021 


0.081 


5.0-5.5 


0.013 


0.062 


5.5-6.0 


0.007 


0.043 


" Redshift. 


b ' c Fractional 


contribution of 



sources in the redshift bin to the total sources 
with 1.1-mm flux density greater than the 
specified limit. 

tribution and re-distributing those galaxies into the isolated disc 
subpopulation. The predicted cumulative number counts are shown 
in Fig. [7] in which the total single-dish counts are also plotted 
for comparison. The values of the counts and the fractional con- 
tributions of the isolated discs and starbursts for various flux bins 
are given in Table [6] To facilitate comparison with observations, 
the table includes approximate flux densities for ALMA Bands 6 
and 7 calculated using the mean ratios for our simulated SMGs, 
Salma-s/Si.i = 0.8 ± 0.06 and Salma-t/Si.i = 1.6 ± 0.2. 

As for the single-dish counts, the isolated discs dominate at 
the lowest fluxes (5*1.1 < 1 mJy, or Ssbo < 2 mJy), and the bright- 
est sources are pre-dominantly starbursts. Because the galaxy-pair 
SMGs are resolved into multiple fainter sources, the cumulative 
number counts for ALMA-detected SMGs are lower at all fluxes by 
~ 30 — 50 per cent, and the differential counts should be steeper. 
(At fluxes fainter than those shown, for which the isolated discs 
completely dominate and the bright sources contribute negligibly 



to the cumulative counts, the single-dish and ALMA cumulative 
co unts will be almost id enti cal.) This effec t has a lso been discussed 
bv lKovacs etail d2010h and lSmolcic et all d2012h . 

The redshift distributions for two flux cuts are shown in Fig. 
[8] and the values are given in Table [7j The mean redshifts for the 
Si.i > 1.5 and 4 mJy (Ssso > 3.5 and 9 mJy) bins are 3.0 and 
3.5, respectively. The mean values are almost identical to those for 
the single-dish counts, and there is also the same tendency for the 
brightest SMGs to be at higher redshifts. The redshift distributions 
are similar, but the distribution for the Si.i > 4 mJy sources is 
less strongly peaked than for the single-dish sources. The latter's 
redshift distribution is more strongly peaked because the redshift 
distribution of the bright galaxy pairs peaks at z ~ 3.5. 



6 DISCUSSION 

6.1 Are modifications to the IMF required to match the 
observed SMG counts? 

One of the primary motivations for this work is to reexamine the 
possi bility that SMG number counts provide ev idence for a flat IMF 
(B05, ISwinbank et alJr2008l; iDave et alj|2010h . To test this claim, 
we have assumed the null hypothesis - that the IMF in SMGs does 
not differ from what is observed locally - and used a Kroupa IMF. 
Furthermore, our model is constrained to match the present-day 
mass function - or, more accurately, to not over-produce z ~ 
massive galaxies - by construction; this is an important 'litmus test' 
for putative models of the SMG population. The counts predicted 
by our model agree well with the observed counts for fields be- 
lieved not to be significa ntly affected by gravitational lensing (e.g., 
lAustermann et alfcOlfjt) . Thus, our model does not require modifi- 
cations to the IMF to match the observed counts. As an additional 
check, we crudely approximate the effect of varying the IMF in 
starbursts in our model by multiplying the (sub)mm fluxes of our 
starbursts by factors of 2 (to represent a 'mildly top-heavy IMF') 
and 5 (to represent a 'flat IMF', because this is the factor appropri- 
ate for the IMF used in B05; G.-L. Granato, private communica- 
tion). 

Fig. [9] shows the range in number counts predicted for this 



© 2012 RAS, MNRAS 000,fTH20l 



SMG number counts, z distributions, & the IMF 15 




S,.i (mJy) 

Figure 9. Predicted cumulative number counts for the unlensed SMG pop- 
ulation as observed with single-dish (sub)mm telescopes, JV(> Si.i), in 
deg — 2 , versus Si.i (mJy), where we have crudely approximated the ef- 
fects of using a 'mildly top-heavy' ('flat') IMF in starbursts by multiplying 
the 1.1 -mm flux of the starburst SMGs by a factor of 2 (5). The region 
filled with cyan lines indicates the range spanned by IMF boost factors be- 
tween 2 and 5. The counts for our standard model, which uses a Kroupa 
IMF. are shown for comparison. The points and grey hatched region are the 
same as in Fig. [4] The starburst and total counts are increased significantly 
by the 'top-heavy IMF', but the isolated disc and galaxy pair counts are 
not affected because their fluxes are not altered. The starbursts dominate 
the SMG population for S\a > 1 mJy (S 8 50 > 2 mJy). For Si. i > 2 
mJy (Sgso > 5 mJy), the 'mildly top-heavy IMF' model over-predicts the 
counts that are thought not to be boosted by lensing, and the 'flat IMF' 
model overpredicts all the observed counts, which suggests that a signifi- 
cantly top-heavy IMF is ruled out. N.B. The lack of a boost in the counts 
for Si.i < 1 mJy is artificial; see the text for details. 

range of IMF variation. This modification causes starbursts to com- 
pletely dominate the counts for Si.i (Ssso) > 1 (2) mJy (as is the 
case for the B05 model), and the predicted counts are significantly 
greater than most of the observed counts for unlensed SMGs(3 (At 
the lower fluxes, the similarity between the starbursts counts for the 
various IMF factors is an artificial effect caused by the limited pa- 
rameter space spanned by our simulation suite. If the suite included 
mergers of lower-mass galaxies, adopting a top-heavy IMF would 
cause them to contribute to the counts at these fluxes and thus boost 
the counts at all fluxes. Furthermore, the strength of the effect at the 
highest fluxes is underestimated because of the aforementioned ar- 
tificial cutoff of our starburst counts at 5*1.1 > 4 mJy; see 15. II for 
details.) Applying the 'flat IMF' factor causes the model to sig- 
nificantly overpredict all the observed counts for Si.i > 2 mJy 
(5850 jSj 5 mJy), and as explained above, the lack of an effect for 
lower fluxes is an artificial consequence of our merger simulation 
suite not including lower-mass galaxies. The clear conclusion is 
that in our model, significant modification to the IMF in starbursts 
is unjustified. 

Admittedly, the above arguments against IMF variation de- 
pend on the details and assumptions of our model. However, there 
is a much simpler argument against the IMF in starbursts being 
flat or significantly top-heavy: as discussed above (and in detail 

11 Amusingly, t he counts predicted for the 'mildly top-heavy IMF' agree 
well with those o flAretxaga et al.l J201 ll) even though no lensing is included 
in our model. However, field-dependent IMF variation seems unlikely, so 
the preferred explanation is still that the lAretxaga et al.l )201 ll) counts are 
boosted by lensing. 



in Hll and H12), there is growing observational evidence that a 
significant fraction of the single-dish-detected SMG population is 
attributable to multiple sources blended into one (sub)mm source. 
In some cases, the SMGs are early-stage merger s wi th components 
separa ted by > 10 kpc (see lEngel et alj|2010l and iRiechers et al] 
12011 alibi for excellent examples). Hydrodynamical simulations sug- 
gest that at such separations, a strong starburst is typically not in- 
duced (HI 2 ); this is consistent w ith observations of local galaxy 
pairs (e.g., IScudder et al .1 12012|) . When the galaxies eventually 
merge, their SFRs will increase by a factor of a few at the least (typ- 
ically significantly more for major mergers, but we will be conser- 
vative) and if the IMF is flat in starbursts, the (sub)mm flux per unit 
SFR produced will increase by a factor of ~ 5. Thus, the (sub)mm 
flux of the galaxies should increase by a factor of > 10 during 
the starburst. The relative duty cycles of the galaxy-pair and star- 
burst phases are similar (to within a factor of a few), so for a given 
mass and mass ratio, the number density of starbursts and mergers 
in the galaxy-pair phase should be similar. Thus, if the (sub)mm 
flux is > 10 times greater in the starburst phase than in the quies- 
cently star-forming galaxy-pair phase, mergers during the galaxy- 
pair phase should contribute negligibly to the bright SMG popu- 
lation, as is the case in Fig. [4] even when a more modest boost (a 
factor of 2) is used. This is in stark contradiction with observational 
evidence. Indeed, the existence of the galaxy-pair subpopulation is 
a natural consequence of starbursts being very inefficient at boost- 
ing the (sub)mm flux of merging galaxies (see the discussion in 
Section [5T2l l. The only means to circumvent this argument is to ar- 
gue that the multiple-component SMGs observed are all starbursts, 
but as explained above, this is unlikely. 

The above argument does not rule out the possibility of a 
mildly top-heavy IMF in starbursts or systema tic global evolution 
of the IMF with redshift (as suggested by, e.g., Ivan Dokkumll2008l : 
iDavel 120081 ; iNaravanan & Pavel 120121) : it only requires that the 
IMFs in starburst and galaxy-pair SMGs at a given redshift be simi- 
lar. Furthermore, the argument does not rule out the possibility of a 
bottom-heavy IMF in starbursts, which has been suggested recently 
for l ocal massive ellipticals, the descendants of high-z starbursts 
(e.g., Ivan Dokkum & Conroy||2010L 1201 ll: IConrov & van Dokkuml 
|20 1 2t iHopkinsluO 1 21) . If the IMF were bottom-heavy in starbursts, 
the scaling between (sub)mm flux and SFR would be weaker than 
what we have shown in Hll and here. The weaker scaling would 
further increase the contribution of galaxy-pair SMGs to the SMG 
population. Unfortunately, the model and observational uncertain- 
ties are sufficiently large that we cannot use the observed galaxy- 
pair fraction to argue for or against a bottom-heavy IMF in star- 
bursts, but the above arguments suggest that a significantly top- 
heavy IMF is unlikely. 

6.2 Differences between our model and other work 

Because we find that, contrary to some p revious suggestions (B05, 
ISwinbank et alj|2008t bave et alj|2010f) . a top-heavy IMF is not 
required to match the observed SMG counts, it is worthwhile to 
examine why our results differ from those works. There are multi- 
ple reasons our results may differ: the cosmological context (abun- 
dances and merger rates), the evolution of SFR and dust mass for 
individual mergers, the radiative transfer calculation, the effects of 
blending, and differences in the observed counts used. We explore 
these in turn now. 
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6.2.1 Radiative transfer calculation 

In HI 1, we demonstrated that the (sub)mm flux density of our sim- 
ulated galaxies can be well parameterized as a power law in SFR 
and dust mass (see Equations 1 1 5landfT6t. If the same relation does 
not hold in the B05 model, then differences in the radiative trans- 
fer may be one cause of the discrepancy between our counts and 
theirs. While we have been unable to compare directly with the B05 
model, we have compared our results with those of a SAM that uses 
a similar radiative transfer treatment ( lBensonll20ll) . We find that 
relations similar to Equations dl5t and dl6t hold for the SMGs in 
the SAM, so it appears that the radiative transfer is not the primary 
cause of the discrepancy even though some aspects of the radiative 
transfer differ significantly [e.g., the geometry used in the SUN- 
RISE calculations is taken directly from the 3-D GADGET-2 simu- 
lations, whereas the geometry assumed by GRASIL is analytically 
specified and azimuthally symmetric; the GRASIL calculations in- 
clude a sub-resolution model for obscuration from the birth clouds 
around young stars, but we opt not to use the corres ponding sub- 
resolu tion model that is implemented in SUNRISE dGroves et alj 
l2008h for the reasons d iscussed in sections 2.2. 1 and 4.6 of HI 1]. 

iDave et all d20 lOh did not perform radiative transfer; instead, 
they assumed that the most rapidly star-forming galaxies in their 
simulation were SMGs. However, as noted previously (Hll), this 
simple ansatz is not necessarily true because of differences in dust 
mass amongst galaxies and the effects of blending. Thus, if dust 
radiative tran sfer were perform ed in or our fitting functions were 
applied to the lDave et alj ( 120100 simulations, the conclusions might 
differ. 



6.2.2 Merger evolution 

Perhaps the time evolution of the SFR or dust mass in the B05 
SAM and our model differs. B05 parameterize the SFR in bursts 
as M* = Mgas.c/r*, where M gas , c is the cold gas mass and r* 
is a SFR time-scale given by t* = max[/dynTd yn , r*burst,min]. 
Here, /d yn = 50, Td yn is the dynamical time of the newly formed 
spheroid, and 7"*burst,mtn = 0.2 Gyr. The major merger shown in 
fig. 1 of HI 1 has M gas ~ 10 1 Mq when the galaxies are at coales- 
cence. Let us suppose that all the gas is cold. Then, the maximum 
SFR possible give the B05 prescription is 10 11 M Q /0.2 Gyr = 500 
Mq yr _1 , ~ 9 times less than that of the simulation. If the dust 
mass is kept constant, Equations U5I > and \16\ imply that a factor of 
9 decrease in SFR results in a factor of ~ 2.5 decrease in (sub)mm 
flux, which would significantly affect the predicted counts. This 
is of course only a crude comparison, but it demonstrates that the 
SFHs of starbursts in the B05 model may disagree with those in 
our simulations, and differences in the dust content may also be 
important. 

Further evidence that the physical modelling of merger- 
induced starbursts may account for some of the discrepancy is pro- 
vided by the different importance of starbursts in the two mod- 
els. In the B05 model, starbursts dominate the submm counts by 
a large margin and contribute significantly to the SFR density of 
the universe; they dominate over quiescently star-forming discs for 
2 > 3. In our model, isolated discs dominate the (sub)mm counts 
at the lowest (sub)mm fluxes and quiescently star-forming galaxy- 
pair SMGs provide a significant contribution to the bright counts. 
Furthermore, in our model, merger-induced starbursts account for 
< 5 per cent of the SFR density of the universe at all redshifts 
( Hop kins et all l2010ch ; this is consistent with the starburst-mode 
contribution inferred from the stellar densities, surface brightness 



profiles, k inematics, and stellar populati on s of observed merger 
remnants dHopkins et al] l2008bl . l2009ij |bl; iHopkins & Hernquistl 
[2013). 

Differences in the SFHs of mergers may also contr i bute t o 
the discrepancy between our model and that of Dave et al. 



It is not clear that the resolution of the IDave et alj j201Gh simu- 
lation (3.75/i _1 kpc comoving) is sufficient to resolve the tidal 
torques that drive merger-induced starbursts. If it is not, the SFRs 
during mergers of their simulated galaxies would be underesti- 
mated, and part of the discrepancy between their simulated SMGs' 
SFRs and those observed could be attributed to resolution rather 
than IMF variation. Furthermore, their wind prescription may artifi- 
cially suppress the SF R enhance ment in merger-ind uced starbursts 
dHopkins etal]|2012h . Fig. 2 of IDave et all feoich shows that all 
their simulated galaxies lie significantly under the observed SFR- 
Mi, relation, so it is perhaps not surprising that they cannot repro- 
duce the SFRs of SMGs, which are a mix of merger-induced star- 
bursts that are outliers from the SFR-M* relation and very massive 
quiescently star-forming galaxies that lie nea r the relation (HI 2; 
iMagnelli et alJl2012l ; lMichalowski etaill 20121) . 

6.2.3 Cosmological context 

The third major component of the models that can disagree is the 
co smologica l conte xt, which includes the SMF and merger rates. 
In iHavwardl we demonstrated that the predicted number 

densities and redshift distribution of isolated disc SMGs are very 
sensitive to the assumed SMF. Thus, it is worthwhile to compare 
the SMF in the SAMs to the observationally derived SMF we have 
used. While a direct comparison of the B05 SMF to those in the 
literature is complicated by th e B05 model's us e of the flat IMF 
in star bursts (see section 4.2 of lLacev et alj|20icl) . lswinbank et alj 
d2008h have shown that the B05 model under-predicts the rest- 
frame A'-band fluxes of SMGs, which suggests the masses of their 
model SMGs are lower than observed. This would be a natural re- 
sult of an under-prediction of the abundance of massive galaxies. 

If the B05 model under-predicts the SMF, then they must com- 
pensate by making the starburst contribution significantly higher; 
they do this by enabling very gas-rich minor mergers to cause 
strong starbursts (in their model, minor mergers, which only pro- 
duce a starburst if the most-massive progenitor has baryonic gas 
fraction greater than 70 per cen t, account for approxi mately three- 
quarters of the SMG population; iGonzalez et al.l201lh and by mod- 
ifying the IMF in starbursts such that for a given SFR, they produce 
significantly greater submm flux. An under-prediction of the abun- 
dance of all massive galaxies and subsequent need to strongly boost 
the starburst contribution would explain why the relative contribu- 
tions of starbursts and quiescently star-forming galaxies differ so 
significantly. 

6.2.4 Blending 

An additio n al rea son for the discrepancy is that neither B05 nor 
IDave et al] d2010h account for the blending of multiple galaxies 
into one (sub)mm source, which can be significant for both merg- 
ing di scs (HI 1, H12) and physically unrelated galaxies dWang et alj 
l201lh . Our model suggest that galaxy-pair SMGs can account for 
~ 30 — 50 per cent of the SMG population attributable to iso- 
lated discs and mergers. The types of sources IWang et al] d201lh 
observed could also be important. Thus, not accounting for blend- 
ing could account for a factor of ~ 2 discrepancy between the pre- 
dicted and observed counts. 
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6.2.5 Revised counts 

Finally, B05 compared with and lDave et al. utilised number 

counts that have since been superseded by surveys covering sig- 
nificantly larger areas. The new counts are as muc h as a factor of 
2 lower than those of, e.g., IChapman et alj J2005h . so the differ- 
ence in counts is another non-trivial factor that can explain part of 
the discrepancy between our conclusion and that of some previous 
works. 

6.3 Uncertainties in and limitations of our model 

Our model has several advantages: the SEM enables us to isolate 
possible discrepancies between our model and other work that orig- 
inate from differences in the dynamical evolution of galaxy mergers 
and the dust radiative transfer calculation rather than more general 
issues, such as an overall under-prediction of the SMF at z ~ 2 — 4. 
By constraining the SMF and gas fractions to match observations 
and including no additional gas in our models, we match the z ~ 
mass function by construction. Because we use 3-D hydrodynami- 
cal simulations combined with 3-D dust radiative transfer calcula- 
tions, we can more accurately calculate the dynamical evolution of 
galaxies and the (sub)mm flux densities than can either SAMs or 
cosmological simulations. Furthermore, we conservatively assume 
a Kroupa IMF rather than invoke ad hoc modifications to the IMF. 

However, our model also has several limitations: first, com- 
putational constraints prevent us from running simulations scaled 
to different redshifts. Instead, we scale all initial disc galaxies to 
z ~ 3. Ideally, we would run simulations with structural parame- 
ters and gas fractions scaled to various redshifts because the vari- 
ation in the galaxies' physical properties with redshift may cause 
the (sub)mm light curves to depend on redshift. In future work, we 
will run a large suite of simulations that will more exhaustively 
sample the relevant parameter space; such a suite could be used to 
more accurately predict the SMG counts and redshift distribution. 
However, because our predicted counts are dominated by galaxies 
at z ~ 3, our conclusions would likely not differ qualitatively. 

Second, to have the resolution necessary to perform accurate 
radiative transfer on a sufficient number of simulated galaxies, we 
must use idealised simulations of isolated disc galaxies and merg- 
ers rather than cosmological simulations. In principle, inc luding 
gas supply from cosmological scales (e.g jKeres et alj2"0 05^) could 
change the results. However, such gas supply is implicitly included 
in our model for the isolated discs and galaxy-pair subpopulations 
because we use observationally derived gas fractions for these sub- 
populations. Because the duty cycles for mergers are calculated di- 
rectly from the merger simulations, which do not include additional 
gas supply, the gas contents, and thus SFRs, of the starburst SMGs 
could be increased if cosmological simulations were used. How- 
ever, because the (sub)mm flux density depends only weakly on 
SFR [a factor of 2 increase in the SFR increases the (sub)mm flux 
density by ~ 30 per cent; HI 1], this would not change our results 
qualitatively. Similarly, differences in the dust content could affect 
the results, but a factor of 2 decrease in the dust mass decreases the 
(sub)mm flux density by < 50 per cent. 

Third, our predictions are affected by uncertainties in the ob- 
servations used to constrain the model. Of particular importance 
are uncertainties in the normalisation and shape of the SMF, which 
are especially significant at z > 4. For z < 3, the merger rates 
predicted by our model are uncertain by a factor of < 3, but the 
uncertainties are higher at higher redshifts. A decrease in the nor- 
malisation of the SMF or merger rates would decrease the counts 



predicted by our model. However, the most interesting conclusions 
would remain: 1. Galaxy pairs would still provide a significant con- 
tribution to the SMG population. 2. Depending on the factor by 
which the abundances are decreased, a mildly top-heavy IMF might 
not be ruled out. However, a flat IMF would still over-predict the 
counts even if the SMF normalisation and merger rates were de- 
creased by > 10 times, which is surely an overestimate of the un- 
certainty. 

Fourth, we do not include the contribution from physically un- 
related (i.e., non-merging) galaxies blended into a single (sub)mm 
source. From observatio ns, it is known that such galaxies contribute 
to the SMG population l lWang etalJboilh . Whether such sources 
contribute significantly to the population is an open question that 
should be addressed by obtaining redshifts for SMGs observed with 
ALMA. If they do contribute significantly, inclusion of this subpop- 
ulation would increase the number counts predicted by our model, 
change the relative contributions of the subpopulations, and po- 
tentially alter the redshift distributions; however, this would only 
strengthen the evidence against a top-heavy IMF. 



Finally, we do not model the effects of gravitational lens- 
ing. Whereas the fainter sources are likely dominated by un-lensed 
SMGs, the bright counts (Si.i > 5 mjy, o r Ssso > 11.5 mjy) 



are likely boosted significantly by lensing i Negrello et al 


|2007, 


l20ld;|Paciga et alJl2009l;lLima et alj|20ld: Vieira et alj|2010 


). Fur- 



thermore, t he likely cause of the dis crep ancy between the cou nts 
observed bv lAustermann et al] d2010h and lAretxaga et all d201ll) is 
that the latter counts are boosted by galax y-galaxy weak lensing 
(see the discussion in lAretxaga et al|201ll). so the d i screpa ncy be- 
tween our predicted counts and the I Aretxaga et"al1 d201 ll) counts 
might be resolved if we included the effects of lensing from a fore- 
ground overdensity of galaxies. Again, our model predictions are 
conservative because inclusion of gravitational lensing would boost 
the counts; thus, inclusion of lensing would strengthen the argu- 
ment against a top-heavy IMF. 



7 CONCLUSIONS 

We have presented a novel method to predict the number counts and 
redshift distribution of SMGs. We combined a simple SEM with 
the results of 3-D hydrodynamical simulations and dust radiative 
transfer to calculate the contributions of isolated discs, galaxy pairs 
(i.e., infall-stage mergers), and late-stage merger-induced starbursts 
to the SMG number counts. Our model is constrained to observa- 
tions as much as possible; consequently, we are able to isolate the 
effects of uncertainties related the dynamical evolution of mergers 
and the dust radiative transfer - which are perhaps uniquely rele- 
vant to the SMG population - from more general issues that affect 
the high-redshift galaxy population as a whole, such as the SMF. 
Furthermore, we have conservatively used a Kroupa - as opposed 
to flat or top-heavy - IMF because we wish to test whether we can 
match the observed counts without modifying the IMF from what 
is observed locally. Our principal results are: 

(i) Our fiducial model predicts cumulative number counts that 
agree very well with those observed for fields thought not to be 
significantly affected by lensing. 

(ii) Except at the lowest fluxes (5*1.1 < 1 mjy, or Ssso < 2 
mjy), merger-induced starbursts account for the bulk of the pop- 
ulation not accounted for by galaxy-pair SMGs, and the brightest 
sources are pre-dominantly merger-induced starbursts. Thus, iso- 
lated discs contribute negligibly to the bright SMG population. The 
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contribution of isolated discs to the SMG population is a robust 
testable prediction of our model. 

(iii) Contrary to the conventional wisdom, bright SMGs are not 
exclusively merger-induced starbursts; our model predicts that qui- 
escently star-forming galaxy-pair SMGs account for ~ 30 — 50 
per cent of SMGs with Si.i > 0.5 mjy (Ssso > 1 mjy). Though 
the precise fraction is sensitive to the details of the modelling, the 
prediction that galaxy pairs contribute significantly to the popula- 
tion (i.e., tens of per cent rather than a few per cent or less) is ro- 
bust. The observational diagnostics presented in H12 can be used 
to determine the subset of the SMG population that are quiescently 
star-forming, thereby testing this prediction of our model. 

(iv) The typical redshifts of the model and observed SMGs are 
similar, but the model may over-predict the number of SMGs at 
2 > 4. This may be because the SMF used in our model over- 
predicts the number of massive galaxies at those redshifts; thus, 
observations of the abundance of SMGs at z > 4 may provide 
useful constraints on the massive end of the SMF at those redshifts. 

(v) Because we have not modified the IMF, our results suggest 
that we cannot reject the null hypothesis that the IMF in high- 
redshift starbursts is no different than the IMF in local galaxies. 
A crude test suggests that if we use even a mildly top-heavy IMF in 
our model, the SMG counts are significantly over-predicted. Thus, 
we conclude that the observed SMG number counts do not provide 
evidence for a significantly top-heavy IMF. 

(vi) There are multiple possible reasons our conclusions differ 
from those of previous work, including differences in the radiative 
transfer calculations, the merger evolution, and the cosmological 
context and the lack of a treatment of blending in previous work. 
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